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ABSTRACT 

Context. We investigate the young (proto)stellar population in NGC 2023 and the L 1630 molecular cloud bordering the H n region 
IC 434, using Spitzer IRAC and MIPS archive data, JCMT SCUBA imaging and spectroscopy as well as targeted BIMA observations 
of one of the Class protostars, NGC 2023 MM 1. 

Aims. We study the distribution of gas, dust and young stars in this region to see where stars are forming, whether the expansion of 
the H II region has triggered star formation, and whether dense cold cores have already formed stars. 

Methods. We have performed photometry of all IRAC and MIPS images, and used color-color diagrams to identify and classify 
all young stars seen within a 22' x 26' field along the boundary between IC434 and L 1630. For some stars, which have sufficient 
optical, IR, and/or sub-millimeter data we have also used the online SED fitting tool for a large 2D archive of axisymmetric radiative 
transfer models to perform more detailed modeling of the observed SEDs. We identify 5 sub-millimeter cores in our 850 and 450 yum 
SCUBA images, two of which have embedded class or I protostars. Observations with BIMA are used to refine the position and 
characteristics of the Class source NGC 2023 MM 1. These observations show that it is embedded in a very cold cloud core, which 
is strongly enhanced in NH2D. 

Results. We find that HD 37903 is the most massive member of a cluster with 20 - 30 PMS stars. We also find smaller groups of 
PMS stars formed from the Horsehead nebula and another elephant trunk structure to the north of the Horsehead. Star formation is 
also occurring in the dark lane seen in IRAC images and in the sub-millimeter continuum. We refine the spectral classification of 
HD 37903 to B2 Ve. We find that the star has a clear IR excess, and therefore it is a young Herbig Be star. 

Conclusions. Our study shows that the expansion of the IC 434 H 11 region has triggered star formation in some of the dense elephant 
trunk structures and compressed gas inside the L 1630 molecular cloud. This pre-shock region is seen as a sub-millimeter ridge in 
which stars have already formed. The cluster associated with NGC 2023 is very young, and has a large fraction of Class I sources. 
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0^ 1 . Introduction (2.2 jum) down to 14 mag. They identified 4 embedded clusters, 

the richest being associated with NGC 2024 with over 300 em- 

- The southern part of the Orion B giant molecular cloud complex, ^^^^^^ ^^^^^ ^^^^^^ associated with NGC 2023, a re- 

L 1630, borders the arge H 11 region IC 434, which is expanding 3^^^^^^^ nebula just north of the Horseh ead with only 21 embed- 

^ , into the molecular cloud and possibly triggering star formation. ^^^^.^^^ ^ ^^^.^ ^^^^^^^ ^^^^ bv lDepov et all ([T990l) in J, 

■ The mterface between the molecular cloud and the H 11 region ^ ^-^^ ^ ^.^.^^ ^ magnitude of 15 red uced the number 

IS seen as a bright north-south ridge of glowmg gas with the ^ ^^^.^ 2023 cluster to 16. iLaunhardt et all 

Horsehead nebula (B33) and several smaller elephant trunks ^^99^ ^^^^^ the CS J = 2 ^ 1 cores for dust emission at 

seen in silhouette against the bright nebulosity. B33 points di- j 3 ^^^^^ ^^^^ associated with compact 

rectly towards t he binary system, cr On, which is ionizing the H n H,ist emission 

region IC434 (lAbergel et al.ll200l iPoundl iReipurth & Ballvl) . ■ 

Since the Horsehead is seen in absorption against the nebula, it In the vicinity of NGC 2023 ILaunhardt et al] 11996') found 

must be at the same distance or cl oser to us than cr Ori, which has faint extended emission associated with the southern part of the 

a distance of 330-385 pc (Bro wn. De Geus & de Zeeuwlll994l: reflection nebula; to the west of the nebula they found a rela- 

ICaballerol l2 5b 8). Although 350 pc is a more probable distance, lively bright double source LBS 34 SM 1 and SM 2 embedded in 

in this paper we adopt the distance 400 pc, which allows easy fainter extended emission. jandeU et al.| (Il999|) serendipitously 

comparison to other studies of this region. The derived luminosi- re-discovered the latter two sources while carrying out a CO 

ties may therefore be overestimated by 20 to 30%, and physical 7 = 2^1 survey of the Horsehead and L 1630 and found 

sizes by 10-20%. an extreme high velocity outflow associated with and driven by 

The first systematic search for dense cores and embedded LBS3 4 SM 1 , which they called NGC 2023 MM 1 . |SandeU et al.] 

star clusters in the L 1630 c loud comp l ex wa s carried out by (119991) caiTied out extensive observations of NGC 2023 MM 1 

iLada. Bally & StarkI (Il991h : iLada et al.l (|1991), who surveyed and MM 2 and found that MM 1 had all the characteristics of a 

3.6 square degrees in CS 7 = 2 ^ 1 (98 GHz) and in K-band Class protostar with a luminosity of < 10 Lq. There are other 

signposts of star forma tion in the vicinity of NGC 2 023 such 
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free emission sources (iReipurth et al.l2004l) without any obvious 
optical or near-IR counterparts. These observations suggest that 
there may be a more embedded population of protostars hidden 
in t he dense molecular c loud cores surrounding NGC 2023. 

iBowler et al.l (l2009l) have recently derived an infrared cen- 
sus of star formation in a 7'x7' region around the Horsehead 
nebula. They used deep near-infrared (IRSF/SIRIUS J H K,) 
and mid-IR (Spitzer IRAC) observations. Their findings support 
triggered star formation due to radiation driven implosion in the 
Horsehead, although they see no evidence for sequential star for- 
mation in the immediate vicinity of the cloud/H ii region inter- 
face. They identify three bona fide and five candidate young stel- 
lar objects (YSOs) in the Horsehead area. 

In this paper we report on millimeter aperture synthesis 
observations of NGC 2023 MM 1 and MM 2 with BIM^H to 
better characterize the protostars and to see if we could re- 
solve its accretion disk. We also carried out more extensive 
continuum mapping of the NGC 2023 field with SCUBA on 
JCM10 to look for young deeply embedded stars, which can 
be too heavily obscured to be detected in the near infrared 
(e.g. NGC 2023 MM 1). We have also retrieved all IRAC and 
MIPS images from the archives of Spitzer Space Observatory; 
these data go much deeper and to longer waveleng ths than those 
availa ble from any ground based observatory. The iBowler et al.l 
( 1200 9) study of the Horsehead area uses the same IRAC data 
as we analyze in this study, but we examine a much larger 
field. Since we started this project, two other S CUBA studies of 
the Orion B South region have been pu blished (iJohnstone et al.l 
120061 iNutter & Ward-Thompsonll2007h . Both these studies are 
in part based on the same SCUBA observations. Although we 
do not go as deep as Johnstone et al. (2006), we have performed 
a very careful data (re-) reduction, to achieve better astrometry, 
sharper images and well calibrated data. This has enabled us to 
identify infrared, radio (and X-ray) counterparts to the sources 
detected in our SCUBA observations. 

2. Observation 

2.1. BIMA observations 

The observations of NGC2023MM 1/MM2 were made on 
Dec 13, 1999 and April 13, 2000 using one frequency setting 
in the C-array configuration of the Berkeley-Illinois-Maryland 
Association (BIMA) interferometric array. The weather condi- 
tions were marginal in December and somewhat better in April, 
but both runs had poor phase stability. The correlator was split 
into four 25 MHz bands giving a velocity resoluti on of ~0.37 
km s~ '. Since the single dish observations by Sandell et alj 
(119991) indicated that NGC 2023 MM 1 has rather strong emis- 
sion in deuterated lines, we decided to look for NH2D 1 1 1 Iqi 
(at 85.9262703 GHz), which is typically quite abundant in cold 
cloud cores. We also included SO 22 — > li, which can probe 
both the outflow and the cold molecular gas in the core. We kept 
the upper sideband free of any molecular transitions that could 
be expected to be excited in a cold low mass cloud core. 3C 273 
was observed for 10 minutes for bandpass calibration and the 
quasar 0609-157 for phase calibration at intervals of 30 minutes. 

' The BIMA array was operated by the Universities of California 
(Berkeley), Illinois, and Maryland with support from the National 
Science Foundation. 

^ James Clerk Maxwell Telescope (JCMT) is operated by the Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. 



The observations of 0609-157 were done using a bandwidth of 
800 MHz for 2.5 minutes. Since there were no planets available 
during the observations we used our phase and bandpass calibra- 
tors for flux calibration. At the time of the observations 3C 273 
had a flux density of 8.8 Jy, while the flux density of 0609-157 
was 5.0 Jy, both with an uncertainty of ~ 10% determined from 
planet observations within a month from our observing runs. 

The data were red uced and imaged in a standard way us- 
ing MIRIAD software dSaultet al.ll995h . The uncertainty in the 
absolute amplitude scale is ~20 %. The data were imaged with 
weights inversely proportional to the variance in order to obtain 
the best signal to noise ratio. We created a continuum map by 
adding together the four windows in the upper sideband, which 
were found to be free of line emission. The synthesized beam for 
the continuum image is 12'.'7 x 7'.'0 at a position angle (pa) - 4°, 
and the continuum image has an RMS noise of 3.6 mJy beam"'. 
We easily detected NGC 2023 MM 1 with a flux density of ~ 
110 mJy. However, MM 1 was found to be more extended than 
what we determined from the JCMT SCUBA images, which 
appears implausible. We therefore also imaged our phase cal- 
ibrator, 0609-157, which should be a true point source in the 
BIMA beam, and found it to have a Full Width Half Maximum 
(FWHM) of 6"x 3", i.e. the beam is broadened by ~ 1'.'5 - 2". 
Adding 2'.'5 in quadrature to the theoretical synthesized beam, 
we find MM 1 to be unresolved. Since 0609-157 was observed 
only every thirty minutes, it is impossible to judge what the true 
synthesized beam size was for our observations. It is very likely 
that if we re-observed MM 1 in perfect sky conditions with the 
same spatial resolution, we would find it at most marginally re- 
solved. In order to check whether the phase errors were so large 
that we would lose photons out of the main beam, we also deter- 
mined the integrated flux over the central I'xl' area centered on 
MM 1 . This gave the same flux density as we obtained from fit- 
ting a Gaussian to MM 1 . We therefore believe that we retrieved 
all the flux, i.e. the error in flux density is dominated by the un- 
certainty of the assumed flux density for 3C 273 and 0609-157, 
both of which are variable. 

We easily detected NH2D In — > loi, but not SO, nor any 
other molecular transition in our data set. The RMS in an indi- 
vidual channel, i.e. for a velocity resolution of 0.34 km s is ~ 
140 mJy beam"'. 

2.2. SCUBA observations 
2.2.1. SCUBA Scan maps 

We obtained several large scan maps (8' x 8') of the NGC 2023 
region wit h the Submillimeter Common User Bolometer Array 
(SCUBA) (i HoUand et al ll 19991) on JCMT during SCUBA com- 
missioning time on four nights in December, 1997. SCUBA has 
37 bolometers in the long (850 fim) and 91 in the short (450 /im) 
wavelength array separated by approximately two beam widths 
and arranged in a hexagonal pattern. Both arrays can be used 
simultaneously. 

For these observations we used the 850 fim and 450 //m fil- 
ter combination in traditional scan-map mode, i.e. with con- 
tinuous scanning in Azimuth with a scan rate of 32" sec"' 
while chopping in Azimuth with a chop throw of 60". This 
mode gives much flatter baselines than the more commonly used 
"Emerson2" scan maps. All four nights (Dec 11, 17, 18 & 20) 
had excellent submillimeter sky conditions with a precipitable 
water vapor level, PWV < 1 mm. We did a pointing check on 
HL Tau or CRL618 before and after each map, and corrected 
for any pointing drift in the post processing of the data. During 
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the first night, Dec 1 1 , we had large pointing drifts (> 4") during 
the scan map, which resulted in smearing of point- like sources in 
the map, and we decided we could not use it. In total we acquired 
5 maps (12 integrations), half of which were scanned roughly in 
right ascension (i.e with the source near transit) and the rest of 
the maps were obtained with a roughly orthogonal scan direc- 
tion, which was achieved by observing the source while it was 
rising. HL Tau was used as our primary calibrator, although cal- 
ibration was additionally checked using Uranus, CRL618 and 
IRCh-10216. The HPBW (Half Power Beam Width) of SCUBA 
in scan map mode with linear regridding to a 1" grid was mea- 
sured to be ~ S'.'O at 450 yum and ~ 14'.'8 at 850 jum from scan 
maps of Uranus. 

The maps were reduced i n a standard way using 
the SURF r eduction package (iJenness & Lightfool Il999t 
ISandell. Jesso p & Jenness 2001). After reducing each 850 fjm 
map separately, the maps were coadded and regridded onto equa- 
torial coordinates. This map was used to determine the position 
for NGC 2023 MM 1 . We then used this position and corrected 
the individual maps for residual pointing errors (shift and add) 
and perfomed a new coadd, where we also weighted each map 
accounting for the rms-noise in each map. The final positional 
accuracy is estimated to be < r.'5 based o n comparison of the po- 
sitions of MM 1 (this paper) and MM 3 dWyrowski et alJl2000h 
deduced from BIMA observations. The rms noise level in the fi- 
nal maps is 45 mJy beam"' for the the 850 yum map, and ~ 0.45 
Jy beam"' for the 450 yum map, 

2.2.2. SCUBA jiggle maps 

We also took several long integration jiggle maps of the 
NGC 2023 MM 1/MM 2 field on Dec 1 1 1997. These maps were 
reduced separately and added to earlier jiggle maps discussed in 
ISandell et all (IT999). A subset of the 850 jum maps (only those 
which had pointing observations immediately before or after the 
map) were used to determine the position of MM 1 . All maps 
were then corrected to this position and coadded by accounting 
for the rms-noise in each individual map. The jiggle maps are 
likely to have lower astrometric accuracy than the scan maps, 
due to the limited number of pointing checks (Table [1]) we did 
for the scan maps. The rms noise level in the jiggle maps is 30 
mJy beam"' at 850 /urn and 0.23 Jy beam"' at 450 /urn. The cal- 
ibration accuracy is estimated to be better than 10% at 850 fim 
and about 15% at 450 fim for both the jiggle and the scan maps. 
The flux densities deduced for MM 1 and MM 2, which were 
observed in both modes, agree exceptionally well, see Table [1] 

All ou r maps were conv erted to FITS-files and exported to 
MIRIAD dSault etal. Ill995h for further analysis. In order to cor- 
rect for the error lobe contribution, especially at 450 jum, we have 
deconvolved all the maps using CLEAN and a circular model 
beam deduced from observations of Uranus. For the scan maps 
we did the deconvolution for only the three sub-fields where we 
have strong dust emission: the MM 1/MM 2 region, NGC 2023 
and MM 3 and MM 4, and the MM 5 region. The results of this 
analysis are discussed later in the paper. 

2.3. Ancillary molecular line data 

We also have access to large '^CO and '^CO J - 2 —> 1 scan 
maps obtained as bad weather backup programme on JCMT be- 
tween 1995 and 1998. The maps were obtained by scanning in 
RA or Dec with sampling every 5" in the scan direction and step- 
ping by 10" in the orthogonal direction. Since the maps were 



obtained in wet weather conditions and over a long time period, 
each sub-map was made to overlap with previous maps to ensure 
that the calibration is uniform over the whole map. Most maps 
were repeated several times, resulting in integration times per 
spectral data points between 5-30 seconds. The '^CO and '^CO 
observations cover an area from south of the Horsehead nebula 
to about the southern part of NGC 2023, and target the L 1630 
cloud boundary towards the IC 434 ionization front. To the east 
these observations stop just short of NGC 2023. The '-CO map 
is ~ ir X 14' (a X 6) and contains ~ 12,000 spectra with a ve- 
locity resolution of 0.10 km s"'. The noise level of the map is 
on the average ~ 0.15 K km s"', although there are areas where 
the noise level can be five times higher. The '^CO map covers a 
slightly smaller area, ~ 9'-10' x 14', has ~ 9,600 spectra and a 
noise level similar to the '^CO map. 

We have also retrieved, reduced, and analyzed the '^CO and 
C'**0 J - 2 -> 1 maps f r om th e JCMT raw data archiveQ 
discussed in iJohnstone et al.l (l2006h . as well as a smaller C'^O 
7 = 2^1 map of flie MM 1/MM 2 cloud core. 

Any comprehensive analysis of these data sets are well be- 
yond the scope of this paper, but the large CO and '^CO maps 
have been very valuable to look for outflow activity or tempera- 
ture enhancements in the vicinity of the sub-millimeter sources. 
Additional details of these data set are therefore given, at those 
places in the paper where we make use of them in this study. 

2.4. Spitzer IRAC and MIPS (24 fim) images 

We extracted observations from the Spitzer Space Observatory 
archive (Program ID 43: an IRAC Survey of the L1630 and 
L1641 (Orion) Molecular Clouds - Fazio et al.). The IRAC 
data were obtained in the 12 sec High Dynamic Range (HDR) 
mode using four different AORs (Astronomical Observation 
Requests), all with overlap. The HDR mode takes pairs of im- 
ages with 0.6 and 12 sec frame times at each position in all four 
bands. Each star in the field is therefore covered with at least 
eight image frames. We have processed both short (0.6 sec) and 
long (12 sec) integration basic calibrated data (BCD) frames in 
each channel using the Artifact mitigation software developed by 
Sean Carey and created mosaics using MOPEX. We used APEX 
to identify and extract photometry of all point sources in both the 
short and long integration images of the 3.6 and 4.5 yum bands, 
while we only used the short integration images for the 5.8 and 
8.0 fj.m band. 

We also extracted MIPS 24 /urn observations for this region 
from the Spitzer Space Observatory archive (Program ID 47: A 
MIPS Survey of the Orion L1641 and L1630 Molecular Clouds 
- Fazio et al.). The MIPS 70 yum image was severely saturated 
and could not be used, but the MIPS 24 fim data were of good 
quality with only two stars in the southern part of NGC 2023 
partially saturated. The MIPS dataset was flatfielded, overlapped 
and mosaiced using MOPEX. 

3. Large scale overview of the NGC2023 

A false-color image of the NGC 2023 region is shown in Fig. [T| 
(blue: 4.5 //m; green: 8.0 //m; red: 24.0 yum). Figure |2] shows 
the 850 yum SCUBA image overlaid as a contour map on the 
Spitzer IRAC 8 fim short integration image. Also marked in 
Fig. |2] are the locations of all stars identified as young stellar 

^ The JCMT Archive project is a collaboration between the Canadian 
Astronomy Data Centre (CADC), Victoria and the James Clerk 
Maxwell Telescope (JCMT), Hilo. 
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In the SCUBA images (Fig. |2|l we detect emission from 
the infrared-dark dust ridge. This dust ridge is visible in the 
more extended (and de eper) 850 fim SCUBA image presented by 
iJohnstone et alJ (l2006h . which shows that the dust ridge extends 
past the Horsehead in the south. (In the IRAC images the dust 
lane is much narrower with very narrow "wiggly" filaments.) 
We also see strong emission from the molecular cloud and PDR 
ridge to the south and south east of NGC 2023 and a fainter dust 
lane connected to the more prominent north-south dust ridge. We 
securely identify five sub-millimeter continuum sources named 
MM 1 to MM 5, detected at both 450 and 850 ^im, all of which 
appear to be in an early stage of evolution. Four of the five 
source s de tected by us were also detected by IJohnsto ne et al.l 
( i2006h and lNutter & Ward-Thompsonl(l2007h using SCUBA ob- 
servations. The sub-millimeter sources MM 1 and MM 2 are lo- 
cated in the northern part of dust ridge. The only sub-millimeter 
source detected in the mid-infrared is MM 3, the bright source 
directly to the south of the NGC 2023. V615 Ori, a heavily red- 
dened T Tauri star lying to the north of the Horsehead nebula, 
though most likely surrounded by an accretion disk, is not de- 
tected in the sub-millimeter However we detect an outflow in 
the '^CO map mostly to the south-east of V615 Ori with rela- 
tively low velocities. 



Fig. 1. False-color Spitzer IRAC and MIPS image of NGC 2023. 
Blue: 4.5 //m; green: 8.0 fim; red: 24.0 //m. The image is centered 
at R.A. = 05''41'"23' and Dec. = -02°19'09", and extends over 
22'x26' (or X J). N is up and E is left. 



objects based on color-color diagrams incorporating 2MASS, 
IRAC and MIPS photometry and the free-free emission sources 
VLA 1-3 identified by Reipurth et al. (2004). The mid-infrared 
and sub-millimeter images provide a broad spectrum overview 
of the stellar content and dust emission from the region. 

The emission in the mid-infrared is dominated by the reflec- 
tion nebula NGC 2023, which is illuminated by the early B-type 
star HD 37903. The images show an infrared dark dust lane at 
a PA of -20°, lying ~ 7' inside the Photon Dominated Region 
(PDR) boundary between IC434 and the L1630 dark cloud. To 
the south-west lies the famous Horsehead nebula, which still can 
be recognized in the mid-IR, although here we only see the sur- 
face layers (PDRs) outlining the nebula. There are a few strong 
mid infrared sources, but overall the images are characterized 
by extended emission from the PDR. This makes identification 
of discrete sources and photometry in the mid-infrared difficult. 
The physical appearance of the region in the mid-infrared dif- 
fers quite significantly from the s tructures seen in low reso- 
lution IRAS and far-infrared data dMookerjea et alj : 



moderate resolution 1.3 mm data (' Launhardt et alj 



20001). the 



1996h and 



the sub-mm SCUBA data presented here. In the far-infrared 
wavelengths the emission is from the HD 37903 and the bright 
sources directly to the south of it. Owing to their poor angu- 
lar resolution, the far-infrared observations detect only a single 
blob of emission midway between the dust emission peak MM 3 
and HD 37903. At still longer wavelengths, dust emission from 
the dark clouds dominates, with only the PDR, illuminated by 
HD 37903, being marginally visible at 850 fim. 



4. Results 

4. 1 . Submillimeter sources in the NGC2023 region 

We have identified a total of five sub-millimeter continuum 
sources in the region mapped with SCUBA. These are the three 
sources MM 1, MM 2 and MM 5, located in the ridge (extended 
along north-south) lying to the west of NGC 2023 and the 
sources MM 3 and MM 4 located directly to the south and the 
south-east of the NGC 2023 nebula. In order to derive the posi- 
tion and sizes of the sub-millimeter sources we fitted a two com- 
ponent elliptical Gaussian using the task IMFIT in MIRIAD, one 
for the sub-millimeter source, and the other for the surrounding 
cloud. The fit to the broader component is mainly to provide a 
good subtraction of the extended emission, and is not to esti- 
mate the flux density of the surrounding cloud. For all the sub- 
millimeter sources the flux densities at 450 and 850 fj.m along 
with positions derived from the 850 fim image and source sizes 
are presented in Table [T] 

Table 1. Results for submillimeter sources mapped with SCUBA 



Source 


^^"2000 


^2000 


5 850 


^450 


e 








mJy 


mJy 


"x" 


MM 1" 


5:41:24.921 


-2:18:08.7 


3.2 ±0.3 


11.8± 1.8 


6.5 X 5.5" 


MM 1 


5:41:24.941 


-2:18:07.9 


3.2 ±0.4 


12.1 ±2.1 


6.7x5.8 


MM 2" 


5:41:26.331 


-2:18:18.6 


1.5 ±0.3 


7.9 ± 1.3 


15.3 x 9.5 


MM 2 


5:41:26.335 


-2:18:17.3 


1.6 ±0.5 


9.6 ±3.6 


19.8 x 11.7 


MM 3 


5:41:37.120 


-2:17:18.1 


1.3 ±0.2 


4.2 ± 0.9 


< 6.6 X 5.6 


MM 4 


5:41:40.520 


-2:17:11.0 


1.6 ±0.4 


11.6±3.4 


30 X 16 


MM 5 


5:41:29.505 


-2:21:18.3 


0.8 ±0.1 


4.5 ± 0.9 


7.1 x5.1 



" Derived from SCUBA jiggle maps. 
Size from 450 ^m 3'.'9 x 3" 



The sub-millimeter emission ridge located to the west of 
the mapped region extends in the north-south direction. There 
is no mid-IR emission associated with this dust ridge. Detailed 
study of the northern part of this ridge using SCUBA im- 
ages at 850 and 450 jim images and CO 7 = 3 — > 2 ob- 
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S''4£»00i' 40" 20' DO* 40* 

RA (J 2000} 



Fig. 2. SCUBA 850 yum contours of the entire mapped region 
near NGC 2023 overlaid on the Spitzer IRAC 8 yum map im- 
age. The peak 850 /im flux is 3.23 Jy/beam and the contours are 
at 0.22, 0.5, 0.75, 1.1, 1.5, 2, 2.5 and 3 Jy/beam. The beam at 
850 yum is shown at the left bottom corner. The crosses indicate 
the location of the sources identified as pre-main sequence ob- 
jects in Fig. [3] as well as the sources presented in Tableland 
the open squares denote the location o f the free-free emission 
sources VLA 1-3 (iReipurth et al.ll2004 . 



servation resulted in the iden t ification of MMl as a class 
protostar (ISandell et all 1 19991) . iJohnstone etaLl (|2006|) also in- 
corporated the SCUBA data presented here together with ad- 
ditional data from the SCUBA archive. The latte r data have 
also been published by 'N utter & Ward-ThompsonI ( l2007h . We 
have completely reprocessed the SCUBA data and have bet- 
ter calibration and astrom etric accuracy cornpared with pre- 
viously published results (t^S andell et al.' '1999^, 'Johnstone et al.' 
l2006l;lNutter & Ward-Thom pson 2007). Launhai'dt et al. (1996) 
detected only the sources MM 1, MM 2 and MM 3 at 1.2 mm 
with the Swedish-ESO Submillimeter Telescope (SEST; beam- 
size of 23") and the IRAM 30 m telescope (beamsize = 12"). 
MM 3, located directly to the south of the nebula is the brightest 
mid-infrared (M IR) source in the region . Surprisingly the coor- 
dinates given by iLaunhardt et al.l (1199 6) for this source are sig- 
nificantly different from those of the 850 and 450 fim sources 
we detect. With the exception of MM 1/ MM 2 all the other sub- 
millim eter s ources wer e detected at 850 irni bv IJohnstone et al] 
( I2OO6I) and iNutter & W ard-Thompson (2007), but their posi- 
tions, sizes and deduced cloud masses diff'er considerably from 
our results. Both of these previous studies identified only one 

su b-millimeter sour ce in the MM 1/ MM 2 cloud core. 

IJohnstone et al.1 (l2006h and iNutter & Ward-ThompsonI 
(|20b7) have identified 14 and 12 sources, respectively, in the 
region we have mapped using SCUBA. We can identify a few 
of their cores in our 850 yum map, but since they are too faint 



to allow secure identification at 450 yum, we do not Ust them 
in Table [1] In Sec. 110.21 we discuss the reasons behind the 
diff'erences between our work and these previous studies. 

4.2. Mid-infrared Ptiotometry using IRAC and MIPS on 
Spitzer 

We carried out multiframe PSF photometry using the SSC- 
developed tool APEX on all the Spitzer IRAC images and on the 
MIPS 24 yum images. While all long integration IRAC images 
are considerably affected by saturation effects, the saturation of 
the 8 yum image is particularly severe and thus could not be cor- 
rected for by the available tools. At 8 yum we therefore used only 
the short integration images to identify sources and derive pho- 
tometry. The nebulous nature of the region and strong emission 
from the associated PDRs makes it extremely difficult to dis- 
entangle the sources from the surrounding clouds and to derive 
accurate photometry. Therefore a combination of automated rou- 
tines and eye-inspection was used to extract the photometry on 
the IRAC and MIPS images. For sources, which APEX failed to 
detect at one or several wavelengths, we used the APEX user list 
option to supply the coordinates for the source to successfully 
derive a PSF fit. 

In all we detected over 1000 stars in the long integration 
3.6 and 4.5 yum IRAC images, while we found only 95 sources 
in the short integration 8.0 yum image. We obtained photome- 
try of 44 sources in the MIPS 24 yum image. The two bright- 
est sta rs, L1630MIR-63 (M M 3), and L1630MIR-73 (stai- # 
219 in IWitt. Schild & Kraim an 1984) in the southern part of 
NGC 2023, are heavily saturated in the MIPS image and their 
flux cannot be measured by conventional photometry. For these 
we used the tool developed by T. Jarrett to rectify saturated stars; 
this tool uses an extended PSF to fit the unsaturated wings of 
the stellar light distribution. In this way we were able to recover 
the integrated flux to an accuracy of ~3-5% based on the re- 
peatability of several measurements from the MIPS BCDs that 
contain the two saturated stars. Since we are mainly interested 
in young stellar objects, we took the 95 sources detected in the 
short integration IRAC 8 yum images as the primary list. We 
cross-correlated this list with sources detected in other IRAC and 
MIPS bands, and also with 2MASS point sources. We used the 
following association radii : 1" for the IRAC images, 2'.'5 for the 
MIPS 24 yum image and 2" for 2MASS data. Table |2| presents 
the coordinates of the 95 MIR sources together with the 2MASS 
magnitudes, Spitzer-IRAC and MIPS flux densities and a pre- 
liminary classification based on selected color-color plots. 

In order to identify the stellar and pre main sequence (PMS) 
stars we used photometry extracted from the IRAC 3.6 and 4.5 
yum long integration images, which go much deeper, have a 
cleaner PSF, and appear to be less affected by nebular emission 
than the 8 yum image. We identify 638 additional sources, de- 
tected in both the 3.6 and 4.5 yum wavebands. Of these sources 
443 are found to have 2MASS point sources associated with 
them. Based on a classification scheme involving the A's-[3.6] 
and [3.6]-[4.5] colors described in Sec. 14.2. 1| we identify 31 of 
the NIR-MIR sources to be Type II PMS stars (Table |3]l. 

We find that, with the exception of L1630MIR-35 and 
L1630MIR-54, afl the MIR sources detected at 8.0 yum were 
detected in the other 3 IRAC bands as well. However only 34 
MIR sources were detected in the MIPS 24 ^m image and 86 
MIR sources were found to have 2MASS counterparts. Out of 
the 44 sources detected at 24 yum, 10 were detected neither in the 
IRAC bands nor in the sub-millimeter images (Table |4|. Six of 
the sources detected only in the MIPS band only are in or near 
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NGC2023. Of the MlPS-only sources (Table H L1630MIPS- 
4 is close to MM2,and L1630MIPS-5 is in the vicinity of 
MM 1 and (probably) the young stars associated with the cold 
MM 1 cloud. Several other MIPS sources lie close to faint 
sub-millimeter peaks that we do not identify as sub-millimeter 
sources in this paper (Fig. |2]l. There are several interesting 
sources (VLA3, Sellgren's star C, HD 37903, and V615 Ori), 
which are bright in the MIR, but which were not detected sub- 
millimeter; see Section|9]for further details. 

4.2.1. YSO classification based on NIR-MIR colors 

We classified the sources detected in the NGC 2023/L 1630 re- 
gion using the NIR and MIR photometric data presented above. 
Figure [3] presents the color-color diagrams derived from the 
2MASS, IRAC and MIPS 24 yum magnitudes of sources detected 
in NGC 2023 along with several criteria (shown as dashed lines 
in Fig.O used to classify them. 

Based on the combined IRACh-MIPS colors ([3.6]-[5.8] vs 
[8]-[24]) (see Figure[3](Ze/f)) we find that almost all the sources 
have significant infrared excess. Most of the sources that have 
been identified in the four IRAC and MIPS bands belong to Class 
II. At 24 jjm the photospheric colors should be close to zero for 
all spectral types, hence the [8]-[24] color is very sensitive to 
excesses. 

For the sources detected in all four IRAC bands we find that 
both the [5.8]-[8.0] and the [3.6]-[4.5] colors for these sources 
extend over 1.5 mag. (Fig. ^middle). Sources with the col- 
ors of stellar photospheres are centered at ([3.6]-[4.5],[5.8]- 
[8.0])=(0,0) and include foreground and background stars as 
well as diskless (Class III) pre-main sequence stars. The box out- 
lined in Fi ^. ^(middle), defines the location of Class II objects in 
the plane (Megeath et al."2004': ' Allen et al.l2004 >). These sources 
have colors that c an be e xplained by you ng, low-mass stars sur- 
rounded by disks. iHartmann et akl (l2005h have shown from their 
observations of Taurus- Auriga young stars that the color criteria 
[3.6]-[4.5]> 0.7 and [4.5]-[5.8]> 0.7 isolate Class O/I protostars 
from Class II and Class III T Tauri stars. 

[Hartmann et al. (2005) have also shown that the color criteria 
Ks-[3.6]> 1.6 and [3.6]-[4.5]> 0.7 discriminate well between 
Class II and Class O/I systems. We have thus used the photometry 
for the sources presented in Table |2] which have 2MASS coun- 
terparts to generate the Kf,-[3.6] versus [3.6]-[4.5] color-color 
plot (Fig.^right). In the same color-color plot we also show the 
31 (out of the 443 sources detected in the 3.6 and 4.5 fim IRAC 
images and with 2MASS associations ) which are identified to 
be PMS in nature based on the color criteria (Table|3]l. We do not 
plot the entire sample of 443 sources in order to avoid crowding 
the figure. Owing to severe contamination due to PAH emission 
at 3.6 yum L1630NIR-23 and L1630NIR-24 do not show colors 
typical of YSOs, but since they are known to be Ha emission 
stars, we count them as YSOs and included them in Table |3] 

A summary of the classification of sources is presented in 
Table |2] For most sources it was possible to arrive at a unique 
class based on all three color-color diagrams. For all sources the 
class derived based on at least two color-color diagrams match. 
For sources indicating two different classes from the three color- 
color diagrams, we have indicated both classes. Out of the 95 
sources detected in the IRAC 8 yum band 32 stars are found to 
be PMS stars (Classes I, II and I/II) based on their colors. Two 
sources detected only in the IRAC 5.8 and 8.0 //m (L1630MIR- 
35 and 54) are most likely PMS, while an additional 14 stars 
marked as type II/III may or may not be PMS stars. The 10 
MlPS-only sources (Table|4]i are likely to be PMS stars, although 



we cannot exclude the possibility of some of these stars being 
background stars. Use of the more sensitive long integration 3.6 
and 4.5 jim images in combination with 2MASS data yields an 
additional 31 PMS objects (Table |3]l. Since the long integration 
3.6 and 4.5 fim images go deeper and are less affected by ex- 
tinction than the 2MASS survey, we also checked how many 
of these IRAC sources are likely to be PMS stars. Since there 
are no strong PAH (polycyclic aromatic hydrocarbon) features 
in the 4.5 fim band, and the 3.3 and 3.4 fim PAH bands will be 
in emission rather than absorption, the [3.6]-[4.5] color should 
be sensitive to stars with a very red continuum, i.e. stars with in- 
frared excess. If we use the color criterion, [3.6]-[4.5] > 0.2, for 
sources detected only in 3.6 and 4.5 fim which have no 2MASS 
counterparts, we find another 38 sources which are likely to be 
PMS stars. In Table |2] we denote the sources which exclusively 
occupy the region characteristic of the stellar photospheres in the 
color-color diagram as "S", although some of them may indeed 
be young stars. 

Figure |2] shows the PMS objects identified in the region with 
vertical white crosses. Thus in the NGC 2023 region we identify 
73 out of a total of 739 sources (1 1 1, if we include the very red 
sources seen only in the IRAC 3.6 and 4.5 ;um images) to be 
PMS in nature and of these five are bona fide Class I sources. 
Nine are identified as Class I/II. We discuss the distribution of 
these YSOs in Sec.fTOl 

Table 4. Sources detected only in the MIPS 24 //m band 



Source 


Q'2000 


^2000 


F24 


LI 630 






mJy 


MIPS-1 


5:40:51.48 


-2:28:48.5 


2.9±0.2 


MIPS-2 


5:40:51.57 


-2:27:49.2 


6.0±0.2 


MIPS-3 


5:40:51.84 


-2:25:53.5 


2.7±0.2 


MIPS-4 


5:41:16.64 


-2:23:58.8 


34.7±0.2 


MIPS-5 


5:41:25.02 


-2:18:15.3 


2.0±0.2 


MIPS-6 


5:41:26.35 


-2:18:21.1 


11.2±0.2 


MIPS-7 


5:41:37.03 


-2:18:17.8 


44.1 ±0.2 


MIPS-8 


5:41:43.75 


-2:16:46.8 


37.4±0.2 


MIPS-9 


5:41:44.36 


-2:17:07.9 


55.9±0.2 


MIPS- 10 


5:41:47.20 


-2:09:40.8 


6.0±0.1 



5. The NGC2023 MMI/MM 2 region 

As mentioned earlier the SCUBA data for the NGC 2023 
MM 1/MM 2 region is of improved quality compared to the data 
presented by S andell et aL (1999 ). We determine the position of 
MM 1 more accurately with BIMA and we resolve MM 2 in both 
SCUBA bands. We do not detect MM 1 with either 2MASS, 
IRAC or MIPS. The flux densities for MM 1 and MM 2 at 450 
and 850 jim (Table [TJ are within 10% of the flux densities mea- 
sur ed bv.Sandell et al..(1999,) . 

ISandell et al.l d 19991) showed that the sub-millimeter source 
MM 1 satisfies all the criteria for a low mass class protostar It 
is not detected in the near-infrared. It powers a jet-like outflow 
with extrem e velocity (Av > 20 km s ') detected in the CO 
7 = 3^2 (ISandell et al.ll 19991) . Its outflow efficiency, defined 
as the ratio of the momentum flux due to the mechanical force 
to the radiative momentum flux (Fco/Z^rad), is > 1000, a value 
typically seen in class sources. A large fraction of its lumi- 
nosity is emitted in the sub-millimeter Lsubmm/i'boi > 10""^ with 
Lboi < 10 Lq. They found the protostellar disk to be extended in 
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Fig. 3. Col or-color diagrams for a ll MIR and PMS-NIR sources. Approximate classification criteria adopted from iMegeath et aTl 
(|2004 and lMuzeroile et alJ (l2004l) are shown respectively for the left and the middle panel, (a) MIR color-color diagram based on 
Spitzer 3-band IRAC and MIPS photometry, (b) 4 band IRAC color-color diagram an d (c) NIR-MIR c olor-color diagram based 
on 2MASS A's and two Spitzer IRAC bands. Dashed lines in all panels are taken from iHartmann et aTl (2005). Reddening vector 
corresponding to the extinction laws for NGC2023/2024 given bv lFlahertvet"al]d2007 ^ is shown in (b). The dashed lines Ks-[3.6] 
= 1.6, [3.6]-[4.5] = 0.7, [4.5]-[5.8] = 0.7 and [5.8]-[8.0] = 0.7 discriminate Class 11 sources from Class 1/0 sources, and ^:s-[3.6] 
= 0.4, [3.6]-[4.5] = 0.2, [4.5]-[5.8] = 0.2 and [5.8]-[8.0] = 0.2 discriminate Class 111 from Class 11. Squares and triangles in (c) 
correspond to sources with and without detection at 8 ;um respectively. 





"2000 "2000 

Fig. 4. Contours of SCUBA 850 and 450 fim continuum jiggle maps of the field surrounding NGC 2023 MM 1 and MM 2 overlaid 
on the grayscale image of Spitzer IRAC 8 pm. The peak intensities at 850 and 450 pm are 3.08 and 9.6 Jy/beam, respectively. The 
contour levels for the 450 pm map are at 0.3, 0.5, 1, 1.75, 2.5, 3.25, 4, 5, 6, 7, 8 and 9 Jy/beam. The contour levels for the 450 pm 
map are at 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 1, 1.5, 2.0, 2.5 and 3.0 Jy/beam. The beam sizes at 450 and 8 50 um are shown at the 
right bottom corner of each panel. "+" and "x" signs mark the sub-mm sources detected in this region by iJohnstone et al] (l2006l) 
and .Nutter & Ward-Thompson (.2007.) . 



the sub-millimeter wit h a size ~ 5", i.e., diameter ~ 2000 AU 
and a mass of 2-4Mo.|Sandell et alj found the dust emission to 
be significantly optically thick at 850 pm. 

5.1. BIMA results 

Figure |5] shows the contours of the 3.5 mm continuum and in- 
tensity of NH2D integrated over the main component observed 
with BIMA, overlaid on the 450 pm SCUBA image. We also 
overplot the contours of integrated intensities of C'^O and C'^O 
7 = 2 — > 1 in Fig. |5] for reference. The sub-millimeter source 
MM 1 is readily detected at 3.5 mm, but not MM 2. NH2D emis- 



sion shows a distribution similar to the continuum emission and 
is detected from both MM 1 and MM 2. Strong emission from 
the deuterated molecule could be a signature of the youth of the 
system. The peak of the NH2D is offset to the north by ~ r.'5 
in declination with respect to the continuum peak at MM 1 in 
the BIMA image. Since the relative positional uncertainty be- 
tween the spectroscopic and continuum observations done with 
BIM A simultaneously is n egligible, the detected offset is signif- 
icant. iRoueff et al.l (l2005l) also found that the NH2D emission 
does not seem to peak at the positions of the embedded proto- 
stars, but instead at offset positions, possibly due to interaction 
between the ambient cloud and the outflow. A plausible explana- 
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Table 3. Coordinates and flux densities of 2MASS & IRAC 3.6/4.5 //m PMS sources in NGC 2023/L 1630. 



Source 
L1630 








•52000 






J 




H 




F3.6 
mJy 


F4.5 
mJy 


Name'' 


NIR-1 


5: 


: 40:46.70 


-2: 


:29: 


:07, 


.1 


16 


.36+0, 


,09 


15.75±0.10 


15.33 +0.16 


0.40 + 0.01 


0.36 + 0.00 




NIR-2 


5 


:40:52.70 


-2: 


:26: 


:25, 


.8 


16, 


.69+0, 


,13 


15.42+0.08 


15.15 +0.14 


0.37 + 0.01 


0.29 + 0.00 




NIR-3 


5 


:40:54.87 


-2: 


:08: 


:33, 


.7 


17, 


.05±0, 


,18 


15.66±0.09 


15.07 +0.13 


0.49 + 0.01 


0.38 + 0.01 




NIR-4 


5 


:40:57.64 


-2: 


:07: 


:13, 


.2 


17, 


.68±0, 


,10 


16.00±0.10 


14.93 +0.14 


0.78 + 0.01 


0.66 + 0.01 




NIR-5 


5 


:40:58.36 


-2: 


:11: 


:34, 


.6 


15, 


.62±0, 


,06 


14.38±0.04 


13.61 +0.05 


2.09 + 0.02 


1.33 + 0.01 




NIR-6 


5 


:40:59.39 


-2: 


:21: 


:24, 


.8 


14, 


.82±0, 


,04 


14.18±0.03 


13.59 +0.05 


2.14 + 0.02 


1.78 + 0.01 




NIR-7 


5 


:41:05.54 


-2: 


:27: 


:17, 


.9 


18, 


.37±0, 


,01 


16.04+0.13 


14.36 +0.07 


1.64 + 0.01 


1.31 + 0.01 


B33-31 


NIR-8 


5 


:41:14.12 


-2: 


:29: 


:01, 


.6 


17, 


.01±0, 


,18 


15.84±0.12 


15.68 +0.26 


0.42+0.01 


0.234+0.01 


B33-33 


NIR-9 


5 


:41:15.61 


-2: 


:14: 


:49, 


.1 


18, 


.70+0, 


,10 


16.21±0.17 


15.10+0.15 


0.44 + 0.01 


0.36 + 0.01 




NIR-10 


5 


:41:18.30 


-2: 


:17: 


:13, 


.6 


16, 


.83±0, 


,15 


15.10±0.06 


14.38 +0.08 


0.72 + 0.01 


0.57 + 0.01 




NIR-11 


5 


:41:19.79 


-2: 


:22: 


:09, 


.2 


18, 


.27±0, 


,10 


16.01±0.13 


14.87 +0.13 


0.50 + 0.01 


0.38 + 0.01 




NIR-1 2 


5 


:41:22.48 


-2: 


:09: 


:15, 


.0 


16 


.93±0, 


,10 


15.74±0.13 


14.60 +0.11 


0.74 + 0.01 


0.69 + 0.01 




NlR-13 


5 


:41:22.49 


-2: 


:17: 


:18, 


.6 


14, 


.68±0, 


,04 


12.95±0.03 


12.05 +0.02 


8.05 + 0.03 


6.89 + 0.03 




NlR-14 


5 


:41:30.06 


-2: 


:14: 


:47, 


.7 


17, 


.67±0, 


,10 


15.54±0.10 


14.42 +0.08 


0.85 + 0.01 


0.76 + 0.01 




NlR-15 


5 


:41:31.93 


-2: 


:17: 


:25, 


.7 


17, 


.19±0, 


,10 


16.03±0.10 


15.00 +0.17 


0.64 + 0.01 


1.18 + 0.01 


HH247 


NlR-16 


5 


:41:32.67 


-2: 


:10: 


:51, 


.9 


16, 


.93±0, 


,14 


15.14±0.06 


14.46 +0.09 


0.72 + 0.01 


0.56 + 0.01 




NlR-17 


5 


:41:32.73 


-2: 


:13: 


:17, 


.5 


17, 


.29±0, 


,23 


15.55±0.10 


14.47 +0.08 


0.82 + 0.01 


0.79 + 0.01 




NlR-18 


5 


:41:32.75 


-2: 


:24: 


:22, 


.0 


18, 


.65±0. 


.10 


16.20±0.15 


14.74 +0.10 


1.01 + 0.01 


0.82 + 0.01 




NlR-19 


5 


:41:32.97 


-2: 


:22: 


:14 


.1 


17, 


.32±0. 


,10 


16.16±0.17 


15.17 +0.14 


0.53 + 0.01 


0.42 + 0.01 




NlR-20 


5 


:41:35.74 


-2: 


:08: 


:55, 


.8 


17, 


.20±0, 


,10 


15.61±0.11 


14.69 +0.10 


0.59 + 0.01 


0.47 + 0.01 




NlR-21 


5 


:41:35.74 


-2: 


:26: 


:35, 


.6 


18, 


.73±0, 


,10 


15.95±0.12 


14.67 +0.10 


0.79 + 0.01 


0.64 + 0.01 




NlR-22 


5 


:41:35.95 


-2: 


:30: 


:16, 


.3 


18, 


.59±0, 


,10 


16.17±0.16 


14.80 +0.11 


0.93 + 0.01 


0.78 + 0.01 




NlR-23'' 


5 


:41:41.51 


-2: 


:15: 


:43, 


.9 


11 


.38±0, 


,11 


11.63+ 0.11 


14.10+0.04 


7.86+0.08 * 


4.00+0.04 


WB 59 


JN1K-Z4 


J 


:41:43.09 


-z: 


:17: 


:39, 


o 
.5 


13 


.03±0. 


,03 


Iz. 35+0.02 


Iz.l 1+0.03 


4.90+0.03 


3.21+0.02 


WH 3 / 


NIR-25" 


5 


:41:44.22 


-2: 


:16: 


:07, 


.5 


13, 


.38±0, 


,04 


12.51+0.05 


11.93+0.04 


14.67 +0.08 


15.03+0.07 


13AV218 


NIR-26 


5 


:41:45.48 


-2: 


:19: 


:46, 


.4 


17, 


.52±0, 


,10 


16.28+0.18 


14.95+0.14 


0.85+0.01 


0.69+0.01 




NlR-27 


5 


:41:46.67 


-2: 


:28: 


:38, 


.1 


16, 


.72±0, 


,14 


15.98+0.12 


15.24+0.17 


0.45+0.01 


0.39+0.00 




NlR-28 


5 


:41:46.70 


-2: 


:17: 


:13, 


.9 


16 


.92±0, 


,16 


16.20+0.18 


15.26+0.18 


0.43+0.01 


0.38+0.01 




NlR-29 


5 


:41:50.44 


-2: 


:17: 


:39, 


.6 


18, 


.20±0, 


,10 


16.39+0.20 


14.56+0.12 


1.11+0.01 


1.01+0.01 




NlR-30 


5 


:41:52.31 


-2: 


:16: 


:11, 


.8 


16, 


.42±0, 


,10 


14.25+0.07 


12.97+0.04 


3.46+0.02 


2.88+0.01 




NlR-31 


5 


:41:56.21 


-2: 


:16: 


:45, 


.9 


16 


.86±0, 


,18 


16.30+0.20 


15.15+0.10 


0.40+0.01 


0.39+0.01 





" Hot emission line star jWeaver & Babcockll2004l) 

* 3.6 fim photometry is unreliable due to severe contamination from PAH emission 

Completely hidden in the Airy pattern of the o f the exceedingly br ight star L1630MIR-72 at both 5.8 and 8.0 jsm. 

Integers c orrespond to NIR sources detect ed bv lDepov et alj | |199(]I) , WB stands for sourc es detected bv lWeaver & Babcocld l l2004h . W are sources 

detected bv lWitt. Schild & KraimanI im4l and B33 are names from lBowler et alj ( |2009|) 



tion could be that NH2D gets destroyed once disks form, so that 
only in the pre-stellar phase the emission is peaked on the pro- 
tostar. For somewhat more evolved protostars it no longer peaks 
at the center of the source. Further, at moderate velocities the 
outflow could sputter NH2D molecules off dust grains, result- 
ing in an enhancement of NH2D. This of course requires that the 
outflow velocities are not high enough to destroy the grains com- 
pletely, but are moderate as found in the shearing layer between 
the outflow and the ambient cloud. 

Although MM 1 appears extended in our BIMA observations 
(Section 12.1b . we do not know whether we resolved the con- 
tinuum emission or not, because of the marginal weather con- 
ditions during our observations. The poor phase stability may 
also affect the position of the continuum emission; q;(2000.0) 
= 05'' 41'" 24'.'929, 5(2000.0) = -02°18' 07'.'2, with an uncer- 
tainty of Ol'S, estimated from dividing up the data into subsets 
and determining the position of MM 1 in each data set. This posi- 
tion agrees within errors with the faint f ree-free emission source 
NGC 2023 VLA 1 (iReipurfli et al.ll2004 . offset by -0'.'0,-Ha.'5 in 
RA and Dec, respectively. Since the free-free emission is only 
0.14 mJy at 3.6 cm, the contribution from free-free emission is 
completely negligible at 3.5 mm, i.e. all the continuum emission 
at 3.5 mm is due to dust. 



5.2. Is NH2D tracing the protostellar disk of MM 1 ? 

Based on CO J=3->2 observations ISandeU et all (Il999h con- 
cluded that NGC 2023 MM 1 drives a strong outflow at a posi- 
tion angle of 145°(measured from north to east). To explore the 
presence or detection of a disk from our NH2D data we assume 
that the disk is oriented at right angles to the outflow with PA 
- 55°. Figure |6] shows the position-velocity (PV) diagram of the 
NH2D emission along the assumed disk as well as the outflow. 
We see a distinct velocity gradient along the disk, suggesting that 
the protostellar disk or the circumstellar material is rotating. The 
emission is slightly redshifted towards northeast and blueshifted 
to the southwest with marginal asymmetry relative to the posi- 
tion of the continuum source which is at an offset of (+ 1'.'2,+ I'.'S) 
relative to the phase center for the BIMA maps. 

Figure |7] shows spectra at selected positions along the disk 
with MM 1 at the center and at the position of MM 2. The spectra 
were analyzed using CLASS method HFS to fit the six hyperfine 
lines. The HFS method of CLASS is a part of the GILDA!^ soft- 
ware package which calculates the total optical depth t, the aver- 
age linewidth of lines and the radiation temperature of the multi- 
plet, assuming LTE. Frequencies and line strengths were adopted 
from Tine et al. (2000) . The derived velocities in the LSR frame. 



* http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 5. The NH2D emission from the MM1/MM2 region is plot- 
ted in grey scale overlaid with contours ofC'^07 - 2 —> 1 
(dotted), C'^O 7 = 2 — > 1 (dash-dot) and 3.5 mm continuum 
(solid black). The extent of the C'^O map, which only covers 
a portion of the image is shown with a grey square. The NH2D 
grey levels are obtained by integrating over the main component 
and are drawn at 0.2, 0.4, 0.6, 0.8, 1, 1.2, 1.4 and 1.6 K km s"'. 
The 3.5 mm continuum emission is centered on MM 1 with con- 
tours drawn at 0.010, 0.025, 0.04, and 0.055 Jy/beam. The BIMA 
beam is shown in the bottom right corner. The and C'^O 
7 = 2 — > 1 emission is integrated over a 1 km s"' wide window 
centered at 9.4 km s"', the systemic velocity of the cloud. The 
contours for C"*0 are from 1% to 100% in steps of 10% of the 
peak value, 7. 1 K km s ' . The contours for C'^O are at 0.7, 0.75, 
0.8, 0.85, 0.9, 0.95, 0.97, 0.99, 1.01, 1.03 and 1.05 K km s'K 
The triangle marks the position of the 3.6 cm free-free emission 
source VLA 1 (coinciding with MMl), while the white cross 
marks the position of the SCUBA continuum source MM 2. 



Vlsr, for the main component, the optical depths and linewidths 
resulting from the fit are presented in Table |5] Based on the fits 
we find that there is distinct shift in Vlsr along the disk, the shift 
is of the same magnitude as the velocity resolution for the BIMA 
observations. Although the velocity resolution of the BIMA ob- 
servations is 0.334 km/s, based on the assumption that the lines 
are Gaussian it is possible to determine the center of the line to 
at least an accuracy of better than a tenth of this value owing to 
the improved accuracy of the fits due to the hyperfine lines. We 
find that the velocity changes by 0.13 km s"' over a lengthscale 
of +3". Most significantly the velocity change per unit length, 
closer to the center of the disk is similar throughout the length 
of the disk. If the velocity gradient was due to Keplerian motion, 
a smaller gradient would have been observed closer to the star. 
This implies that the rotating disk which probably lies hidden 
inside is not seen in NH2D, the emission appears to be primarily 
due to the large scale infalling envelope rotating approximately 
as a rigid body. 

The narrow line widths seen in NH2D also support this non- 
detection. If NH2D were probing a rotating disk, the emission 
from the disk in our relatively broad beam, should by itself 
broaden the line. NH2D is quickly destroyed in shocks and in 



warm gas, so that it may be absent in the accretion disk and what 
we see is the pristine cold gas in the surrounding, slowly rotating 
infalling envelope. We note that although NH2D could exist in 
the dense cold mid-plane of the disk, the size of the disk and the 
low temperatures of the mid-plane would make detection of the 
emission from this region extremely difficult. 

We find that in comparison to the direction along the disk 
there is hardly any velocity gradient along the direction of the 
outflow (Fig.|7]i. This implies that the spread in velocity seen in 
the P-V diagram (Fig. |6]l is primarily due to the broadening of 
the spectral lines by the outflow and the difference in velocities 
of MM 1 and MM 2 (Tab.|5]l. 





10 20 
Velocity (km/s) 



Velocity (km/s) 



Fig. 7. NH2D specti-a along the disk (Left: PA = 55°) at MM 1 
and the spectrum at the position of MM 2 (Right). The green 
curve in each panel shows the results of the hyperfine structure 
fitting obtained using the HPS method in GILDAS/CLASS. In 
each panel the velocity of the main component of the hyperfine 
lines as well as the fitted Gaussian widths are given. 



5.3. Mass of MM 1 and MM 2 from greybody SED fits 

Figure[8]shows the greybody models fitted to the observed SEDs 
of NGC 2023-MM 1 and MM 2. Based on the observed sizes (6,) 
of the sources, for the greybody fitting we have assumed 6^ to be 
5" and 10" for MM 1 and MM 2 respectively. The dust tempera- 
tures (Td), dust emissivity indices (J3) and masses calculated for 
the best fit model are presented in Fig. [8] The total dust and mass 
gas, Mtot has been derived from the fitted and /? using 



,/1200\^^^ 
M,<„ = 1.88X 10"- — — Syie" 



(1) 
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Fig. 6. Position velocity diagram of the NH2D emission along (a) the assumed disk orientation (PA=55°) and (b) along the outflow 
(PA= 145°) centered at MM 1 .The vertical lines denote the six hyperfine spectral lines with yLSR=9.4 km s ' for the main component 
derived from Gaussian profile fitting. 
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Fig. 8. Single temperature greybody fits to the SEDs of 
NGC 2023 MM 1 (left) and MM 2. The fitted dust temperature, 
I'd, the dust emissivity index (fi) and the total mass, Mtot for the 
best fit model are given. 



where d is the distance in kpc, Sy is the total flux (in Jy) at 
frequency v and Mtot is given in [Mq]. In this equation we have 
adopted a gas-to-dust ratio of 100 and standard "Hildebrand" 
mass opacit y, Kq, defined at 25 //m (1200 GHz), i.e. /cnooGHz - 
0.1 cm^g-i (IHildebrandll983h . 

The /?-index de rived for MM 1, 0.7, is much smaller than 
what was found by ISandell et al.l (Il999l) . However, since their 
fit was dominated by single beam sub-millimeter photometry, 
which includes emission from the surrounding cloud, their y6- 
index was too high, and their fit underestimated the SED at long 
wavelengths. We subtract out the emission from the surrounding 
cloud, and since we now also have a flux density measurement 
at 3.5 mm, we obtain a very good fit. A j3 of 0.7, though ap- 
parently low for a Class protostar ( lFroebrichll2005h . is not un- 
usual. HH24 MMS, another Class sour ce in the L 1630 dark 
cloud complex, h as a very similar yS-index (IChandler et al.ll2005t 
iKang et al.ll2008b . We estimate flie total mass of MM 1 to b e 
1.04 Mq, lower than the mass estimate of ISandell et al.1 (1 19991) . 
primarily because we derive a much lower dust emissivity. 

Given that MM 1 drives an energetic outflow the de- 
rived seems quite realistic, while the value of yS derived is quite 



low. As discussed bv lMookeriea et akl (|2007|) both theoretically 
and observationally there is a large scatter in the values of 
yS. Based on current observations a low value of yS typically 
corresponds to evolution of dust grains in circumstellar disks 
around low mass protostars. Since the evolution of grains can 
change their shape, size and chemical composition, it is diffi- 
cult to exactly say wha t causes the lower values o f observed 
yS (B eckwiflietal.H2000h . iMivake & Nakagawal ([199 3) had the- 
oretically estimated that grainsizes of ~ 1 mm correspond to 
yS < 1 . Models of grain growth by collisional coagulation pre- 
dict that while grain growth could occur on very short timescales 
within 10^ years these grains are removed very efficiently within 
a short timescale. This however is inconsistent with observed 
SEDs of classical T Tauri stars. Hence models also suggest that 
small grains probably get replenished, via aggregate fragmenta- 
tion through high-sp eed collision (Dullemond & Dominik 2QQ^ 
iBimstiel et al.l 120091) . We propose that the cause for the mea- 
sured low /3-index is a partially optically thick inner disk, which 
would mimic the appearance of a low /? value. This is plausi- 
ble since many T Tauri sta rs have op tically thick inner disks 
(iBeckwith et al.|[T990: Andrews & WilUams 2005. 2007,) . Hence 
the young disk around MM 1, with is far more massive than 
a Class II disk, could well be partially optically thick at mm- 
wavelengths. 

We note that NGC 2023 MM 2 has a much steeper SED than 
MM 1. Our isothermal fit gives a dust emissivity, /? = 2.08, and 
gives Td = 1 1 .4 K, suggesting that it is a cold pre-stellar core. It 
has a mass of 7.4 M©. Although there is a faint 24 /im source, 
L1630MIPS-4 (Table H ~ 5" from the peak of MM 2, there is 
no sign of activity in the MM 2 core. The 24 jim source is most 
likely unrelated to MM 2, although it could be a low-mass pre- 
main-sequence object embedded in the MM 1/MM2 cloud core. 

We have estimated the mass of the cloud harboring MM 1 
and MM 2 using both the jiggle and scan maps observed with 
JCMT. For this purpose we have smoothed the maps to a reso- 
lution of 30" and then fitted a single Gaussian to estimate the 
total emission from the cloud and MM 1 and MM 2 at 450 and 
850 fiva. The flux densities of the cloud from the scan maps and 
the jiggle maps were determined by subtracting the flux densities 
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of MM 1 and MM 2 in the two types of maps (Tab.[TJ. Finally the 
450 and 850 yum flux densities of only the cloud were obtained by 
averaging the flux densities obtained from the scan and the jiggle 
maps. From the SCUBA observations w e estimate S450 and Sgso 
to be 58.2 Jy and 7.25 Jy respectively. iLaunhardt et all (1 19961) 
using 1.3 mm observations at a resolution of 30" obtained S1300 
to be 1.87 Jy. Using single temperature greybody models we de- 
rive a temperature of 16.5 K and an emissivity index fi - 2.0. 
Using these parameters we estimate the mass of the cloud in 
which MM 1 and MM 2 are embedded to be ~ 18 Mq. 

Thus, from the present SCUBA observations we derive 
a total mass of ~ 26 M0 for th e entire cloud includ - 
ing MM 1 and MM 2. In con trast, iJohnstone et alj ( l2006h : 
iNutter & Ward-ThompsonI (|2007|) . who both identified only one 
core in the MM 1/MM2 cloud complex, using SCUBA studies, 
derived a total mass of 14.7 and 12.7 Mq, respectively, for the 
entire region. 

5.4. Depletion and Deuteration scenario in MM 1 and MM 2 

In order to derive some of the relevant physical parameters, we 
analyzed the NH2D spectral profiles using the HFS method in 
CLASS. This method derives the total optical depth t and aver- 
age width of the lines in an LTE approximation alon g with the 
radiation temperature of the multiplet, Tr. Following lTine et alj 
(1200 0) we have estimated the excitation temperature T^x for the 
multiplet using the following equation 



Tr = [7v(rex) - /v(rbg)](i - e-') 



(2) 



where Jv(T) is the radiation temperature of a blackbody at 
temperature T and Tbg = 2.7 K (the temperature of the cosmic 
background radiation). 

The total column density for NH2D is derived using the op- 
tical depth, excitation temperature and average FWHM Av as 
follows 



^^yll Q(Te.) 



Av- 



SuAul gkTa — 1 



Tul 



(3) 



where v„; is the frequency of the unsplit transition. The val- 
ues of the partition function, Q(Tex), used for the column den- 
sity calculations, ar e 4.47 and 9.37 at 5 K and 10 K respectively 
jRoueffet al.^llool). 

Table |5] also presents the estimated Tex and the column den- 
sities calculated for Tex = 5 K and 10 K. The method HFS as- 
sumes a beam filling factor of unity, this possibly explains the 
rather low excitation temperatures that we derive at all positions. 
The NGC 2023 MM 1 and MM 2 cloud core is extremely rich in 
deuterium, since we observe optically thick emission at many 
positions (Table |5]l. The optical depth and the column density of 
NH2D along the disk shows a centrally peaked profile, with the 
peak being offset from MM 1, while along the outflow the NH2D 
column density increases towards the southeast. However keep- 
ing in mind the assumptions and uncertainties in the estimation 
of these column densities variation by factors of 2 to 3 are most 
probably not significant. Overall, the column densities are sim- 
ilar in ma gnitude ( ~ 10'"* cm"^) to the numbers found in dark 
cores bv iRoueff et al.1 (12005). 

The observed high abundance of singly deuterated ammonia 
in MM 1 is suggestive of freezeout of molecules like CO, since 
depletion of the heavy molecules results in an increase in the 
[H2D+]/[H+] and molecular D/H ratios. Using C"*0 and C'^O 
7 = 2 — > 1 data from the SCUBA archives we have investi- 
gated the depletion scenario in MM 1 and MM 2. The C'**0 data 



(dashed contours in Fig. |5]) does not show any enhancement of 
emission at the positions of MM 1 and MM 2. The C'^O emis- 
sion appears to peak at a position in between MM 1 and MM 2. 

We estimate the molecular hydrogen column densities, 
N(H2), for both MM 1 and MM 2 using the C'**0 and C'^O ob- 
servations. For this we have assumed that the UTE approxima- 
tion holds and have assumed temperatures of 20 K and 1 1 .5 K 
for MM 1 and MM 2 respectively. We note that the estimated 
column density changes only by a factor of 2 for a change in 
temperature from 10 to 100 K. The choice of gas temperatures 
is based on the dust temperatures estimated from the greybody 
fits (Sec. El- Using C'^O (C'^O) we estimate N(H2) to be 
1.6 10^2 cm-2 (2.8 10^2 cm^^) in a 22" beam for MMl. For 
MM 2 we estimate the N(H2) to be 1.7 10^^ cm^^ (3.0 10^2 cm'^) 
from C"*0 (C'^O.) The discrepancy in the N(H2) calculated 
from C"*0 and C'^O column densities is due to C'^O being op- 
tically thick. The relative abundances of NH2D at the positions 
of MM 1 and MM 2 are ~ 10"'' (Table©. 

Alternatively N(H2) can also be derived from the 450 and 
850 //m fluxes observed at MM 1 and MM 2. Using the observed 
flux densities and dust temperatures from the greybody fits we 
get N(H2) in the 22" beam to be 1.2 lO^^ cm"^ 2.6 lO^^ cm"^ 
for MM 1 and MM 2. From dust emission within a beam of 9'.'4, 
comparable to the beam of NH2D, we estimate N(H2) to be 
2.3 10^^ cmr^ and 3.6 lO^^ cm^^ ^j. y^y^ i MU2 respec- 
tively. Thus N(H2) estimated from dust continuum emission is 
larger by a factor of 10 than the N(H2) estimated from the CO 
isotopomers. This is a clear evidence of depletion of CO and its 
isotopomers by a factor of at least 8-10. This is similar to the 
depletion factors observed in prestellar cores, which lie between 
5.5 and 15.5 (.Bacmann et al...200Z) . 

6. The molecular cloud south of NGC 2023 

The reflection nebula NGC 2023 is bounded to the south and 
southeast by a dense, extended molecular cloud coinciding with 
the ridge of strong PAH emission seen in the IRAC and MIPS 
images. At 850 yum the cloud has a size of ~ 98" x 78" PA ~ 
34° measured at the 3-cr contour level with an integrated flux of 
~ 16.7 Jy. Some faint, diffuse dust emission is also seen toward 
the reflection nebula. The strongest dust emission is just out- 
side or overlapping with the bright PDR ridge seen in IRAC and 
MIPS images and the ionize d carbon (C^) ridge. The latter was 
mapped by Wvrowski et all (2000) with the VLA in the C91a 
recombination line. At high resolution t he dust emission breaks 
up i nto several cores (F igure [TOli. Both IJohnstone et"an (l2006h 
and i Nutter & Ward-Th ompson (2007) identified four cores in 
this cloud. We can securely identify only two at both 450 and 
850 jum. One of them (MM 3) is compact and coincides with a 
bright mid-IR source, which has the colors of a Class I/O source 
(Section 14.21 ). The other core, MM 4, is much more extended 
and has no near- or mid-IR counterpart. Even though we see 850 
/vm emis sion at the appro ximate positi ons of the other two cores 
seen b y iJohnstone et alj (12006) and jNutter & Ward-ThompsonI 
( l2007h . they have no counterp arts at 450 /im, and we d o not con- 
sider them as cores. Since the IJohnstone et all (l2006 l) map goes 
deeper, our failure to see these cores could be due to our li mited 
sensitivity. It is, however, p ossible that what lJohnstone et alj and 
INutter & Ward-ThompsonI identify as cores, are structures that 
break up into smaller clumps at the higher angular resolution at 
450 //m, especially since neither core can be identified in any 
molecular line map. 

The molecular gas in the P DR region is hot. The '^CO 
7 = 2 — > 1 map shown in the IJohnstone et alj (l2006l) paper 
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Table 5. CLASS HFS fitting of selected NH2D spectra. MM 1 is at offset (1'.'2,1'.'8). 
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and re-reduced by us from the JCMT archive show '^CO bright- 
ness temperatures of 50 K towards MM 4, which means that the 
gas temperature is > 50 K. Such high gas temperature agrees 
well with the color temperatures of 40-5 K derived fro m far 
infrared observations at 50 and 100 //m (H arvev etalJ[T980 ) one 
arcminute south and west of HD 37903. The color temperatures 
determined between 160 and 100 jum are some what lower, 30-40 
K. The latter color temperatures agree well with lMookeriea et"an 
(I2OOOI) . who determined dust temperatures of ~30 K from 205 
and 138 /im over most of the cloud south of NGC 2023. The hot 
gas is probably restricted to a narrow zone between the reflec- 
tion nebula and the molecular cloud. If we assume a temperature 
of 30 K for the cloud as a whole, we find a cloud mass of 30 
Mq from C'^O J = 2 — > 1, while the mass seen in 850 yim dust 
emission, which only probes the densest parts of the cloud core, 
is only ~ 14.3 M©. Both mass estimates are for the continuum 
and line fluxes from a region extending over 98"x 78". The total 
continuum emission over this region is 16.7 Jy, from which 1.3 
and 1.6 Jy due to the emission from MM 3 and MM 4 are sub- 
tracted to get a flux of 13.8 Jy for the cloud only. The average 
dust temperature is assumed to be 30 K, i.e. similar to what we 
assumed for C'^O. In this case the cloud mass derived from C'^O 
is about twice as large than what we derive from dust emission. 
If C"*0 is undepleted, we would expect to derive a larger mass 
from the C'^O emission than from dust emission, since C'^O 
is also excited in low density gas, while the dust emission only 
traces the very dense gas in the cloud core. If the gas tempera- 
ture is even higher, which could well be the case, it would make 
the difference in mass even larger This strongly suggests that in 
contrast to the MM 1/MM2 region, C'^O is not depleted in the 
MM 3/MM 4 region of NGC 2023. 

7. The NGC2023 MM 3/MM 4 region 

7A. MM 3 

MM 3, coi ncides with an extre mely red IR source, originally dis- 
covered bv lDepov et al.l(ll990h . #3 in their list of near-IR sources 
associated with NGC 2023. The source is also in the 2MASS 
(although not listed by SIMBAD) and the MSX point source 
catalogues. W e detecte d it in all IRAC bands and in the MIPS 
24 ij.m image. IWyrow ski et al. (2000) detected MM 3 at 3 mm 
with a flux density of 36+7 mJy, although they interpreted the 
dust emission to originate in an extended cloud core surround- 



ing the near-IR source detected bv lDepov et al.l (Il990h . We find 
the dust emission to be more compact; at 450 yum the emission 
appears unresolved. The 850 yum size quoted in Table[T]is there- 
fore an upper limit. The excellent spatial agreement between our 
sub-millimeter position for MM 3 and the 2MASS and IRAC 
position (< 0'.'5) therefore supports the classification from the 
IRAC colors, which suggest that MM 3 is a Class I/O protostar. 

E ven though MM 3 is not associated with free-free emis- 
sion dReipurth et al.l l2004t) . it may excite the HH objects 
HH 5, HH4, and HH 1 A- C, which were discovered by 
Mali n. Qgura & WalshI (|1987|) . Examination of the CO 7 = 
2 — > 1 map shows faint high-velocity wings near the position 
of MM 3 and faint red-shifted wings to the east of the source. If 
the outflow is bipolar the blue-shifted molecular outflow would 
therefore be to the west, i.e. in the direction towards the HH ob- 
jects. There is, however, considerable velocity structure in '^CO 
in the NGC 2023 PDR region, and this is by no means a clear 
identification. L1630MIR-62 (Sellgren's star C, see Section m 
is somewhat better aligned with the HH objects. Another pos- 
sibility is a very red near IR source, L1630NIR-15 , which is 
within 1" ofHHl A. 

Because MM 3 lies within the boundary of the reflection 
nebula, which is very bright in the far infrared, there are no 
flux density estimates or even stringent upper limits at 100 
fim, which could help us constrain the temperature of the dust 
emission. However, MM 3 is a visible IR object, presumably 
warmer than 20 K (MM 1) but probably not too hot. We 
therefore assumed a dust temperature of 25 K and estimated 
the dust emissivity and dust mass with a single temperature 
greybody fit using the observed sub-millimeter and millimeter 
flux densities. We obtain a dust emissivity index of 0.8 and a 
dust mass of 0.33 Mq. Although the dust emissivity is on the 
low side, most protostellar objects show dust emissivities in 
the range of 0.7-1.5, and often similar "disk" masses as wel l 
jPent. Matthews and Ward-ThompsonI 119981: iFroebrichI |2005|) . 
re-inforcing our hypothesis that the dust emission originates in a 
disk/envelope around a young protostar. 

Below we analyze in more detail what we can learn about 
MM 3 with the assumption that the dust emission originates in a 
disk or compact envelope surrounding the protostar 
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Fig. 9. SED fits for NGC 2023 MM 3, HD 37903 and V615 Ori. The filled circles indicate the measured fluxes and uncertainties. For 
MM 3, the thin black line represents the best-fitting SED, and the gray lines represent the all other acceptable (x^~Xhesi ^ YSO 
fits. For V615 Ori the thin black line shows the best fitting model which reproduces the spectral class of the source. The dashed 
line shows the stellar photosphere corresponding to the central source of the best fitting model, as it would look in the absence of 
circumstellar dust (but including interstellar extinction). For HD 37903 the model for a stellar photosphere corresponding to a B2V 
star is shown to conclusively demonstrate the infrared excess beyond 4 yum. 




Fig. 10. Contours of SCUBA 850 and 450 yum continuum maps of the field surrounding NGC 2023 MM 3 and MM 4 overlaid on the 
grayscale image of Spitzer IRAC 8 jum. The peak intensities at 850 and 450 yum are 1 .63 and 6.47 Jy/beam respectively. The contour 
for the 850 jum map ai-e at 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 1, 1.2, 1.4 Jy/beam. The contours for the 450 fj.m map are at 1.2, 1.75, 
2.5, 3.25, 4, 4.75, 5.5, and 6.25 Jy/beam. The beamsizes at 450 and 850 ^m are shown at the right bottom corner of each pan el. "+ " 
and "x" signs mark the sub-mm sources detected in this region bv .Johnstone et al.. (.2006.) and Nutter & Ward-Thompson (i2007l) . 
respectively. 



7.2. Nature of MM 3 from SED fitting 

For a more detailed model which explains the emission from 
MM 3 over a wide range of wavelengths we have com- 
bined the SCUBA data with the 2MASS, Spitzer (IRAC, 
MIPS, a nd IRS) and MSX da ta and measurements at 3 mm 
(BIMA). lRobitaiUe et al] (l2007h have created an archive of two- 
dimensional (2D) axisymmetric radiative transfer models of pro- 
tostars calculated for a large range of protostellar masses, ac- 
cretion rates, disk masses and disk orientations. This archive 
also provides a linear regression tool which can select all model 
SEDs that fit the observed SED better than a specified x^- 
Each SED is characterized by a set of model parameters, such 
as stellar mass, temperature, and age, envelope accretion rate, 
disk mass, and envelope inner radius. We have used this on- 
line tool to generate models which fit the observed SED for 
NGC 2023 MM 3. The models assume that stars of all masses 
form via free-fall rotational collapse to a disk and accretion 



through the disk to the star; thus the envelope accretion rate sets 
the envelope mass. Since the distance to the source is somewhat 
uncertain, we have explored a range of distance between 0.3 and 
0.5 kpc and Ay between and 100 mag. Aperture radii used for 
the fits are 30" for the MSX data, 5" for the IRAC and IRS data, 
9" for the MIPS 24 yum data, 14" for the SCUBA data and 13" 
for the 3 mm BIMA data. Left panel of Fig.|9]shows the SED of 
the best-fitting model for MM 3 and one more model which fits 
the data withx^-xl^.t < 10. 

We have obtained the average values and variances of the 
fitted parameters by taking a weighted average of the two mod- 
els for which < 10, in a manner similar to lSimon et al.l 

(,2007) . Both these models correspond to a distance D = 350 pc. 
Here we present the parameters obtained from the two best- 
fitting models, viz., the stellar mass (M»), temperature (r») 
and luminosity (L,), the envelope mass Menv, disk mass Mdisk, 
and the envelope accretion rate Menv The numerical values are 
: M, = 0.88 + 0.22 Mq, L.=21 + 1.4 Lq, r,=3930±168 K, 
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Menv = 0.24+0.02 Mo, Mdisk=0.007+0.004 Mq, and Me„v = 
3.3 10"^ + 8.0 10"^ Mq yr"'. The envelope mass matches rea- 
sonably well with the mass estimated from the simple greybody 
fit. The SED fit is very well constrained in the MIR wavelengths 
owing to the IRS datapoints. As a result we are able to arrive 
at only two models (as opposed to a family of models), which 
reproduce most of the observed features reasonably well, with 
the best fit model fitting the longer wavelength data significantly 
better We note that there is a substantial discrepancy in the flux 
densities measured in the MSX E-band and MIPS 24 /im, al- 
though both bands have very similar wavelengths. The best-fit 
model presented here underestimates the MIPS 24 /im flux. We 
have checked that the quality of the fit does not change signifi- 
cantly by not considering the MSX fluxes, since in the MIR the 
fit is primarily constrained by the IRS fluxes. Thus, only more 
datapoints in the far-infrared region will enable us to arrive at 
a better constrained fit particularly at the longer wavelengths of 
MIR to FIR. 

7.3. MM 4 

Based on '^COJ= 2^1 brightness temperatures > 50 K, MM4 
appears warmer than MM 3. However in the absence of far- 
infrared fluxes it is not possible to derive the dust temperature 
by performing an isothermal fit to the 450 and 850 fim fluxes 
unless the value of /3 is accurately known. We therefore con- 
sidered two isothermal fits constraining to 30 K and then 
to 50 K. In both cases we obtained a fitted value of 1 .5 for 
The masses corresponding to temperatures of 30 and 50 K are 
0.87 Mq and 0.46 M© respectively. Irrespective of the temper- 
ature, the value of j3 for MM 4 is lower than /3 - 2, which is 
typically seen in colder pre-stellar cores. This could be due to 
heating of the core from the outside, which would suggest that 
the core is not well described by a single temperature. Another 
possibility is that the dust has undergone processing and is dif- 
ferent in the PDR region. Which explanation is the right one, 
or whether it is a combination of both, will have to wait until 
we can observe MM 4 with higher spatial resolution at several 
wavelengths in sub-millimeter and far infrared. 

8. NGC 2023 MM 5 region 

MM 5 is a compact object in the southern part of an arc or 
a shell-like structure in the southern part of the dust ridge 
(Fig. fTTT i. There are no optical, infrared or radio sources asso- 
ciated with MM 5 region. Examination of our '-CO and '^CO 
7 = 2 — > 1 maps show no distinctive features in this area 
and the '^CO brightness temperature is ~ 15 K, suggesting that 
the dust emission is associated with a cold prestellar core. An 
isothermal greybody fit to the integrated flux densities suggest 
a dust temperature of 12.2 K, /3 - 2.2, and a mass of 4.3 M©. 
This core was also picked up b y iJohnstone et al.l ( |2006|) and 
iNutter & Ward-ThompsonI (12007 V who find approximately the 
same mass. Both identified the dust emission in the northern 
part of the arc as a separate cloud core (Fig. [TTT l. IJohnstone et al.l 
(1200 6') identified an additional three sub-millimeter "clumps" in 
this region. It is impossible to judge whether all of these are true 
cores, without observing the region in a molecular line tracer 
sensitive to dense, cold gas. 

9. Other sources of interest in the NGC 2023 region 

Comparison of high resolution mid-infrared and sub-millimeter 
data with medium resolution far-infrared data suggests that 



the emission from LI 630, the molecular cloud associated with 
NGC 2023 is primarily extended and diffuse emission. However 
we also find a number of embedded Class I and II objects (see 
Sec. |4.2t , some of which we would have expected to detect in the 
sub-millimeter, since they have clear IR excess and at least some 
are expected to be surrounded by a circumstellar disk and/or en- 
velope. Yet we detect only one MIR source in the sub-millimeter 
Below we highlight some of the stars for which our IRAC and 
MIPS photometry provide new information about their nature, 
viz., HD 37903, V615 Ori, Seflgren's star C, and NGC 2023- 
VLA3. 

HD 37903 or L1630MIR-66 illuminates the reflection neb- 
ula NGC 2023 located in the northern part of the L 1630 molec- 
ular cloud. SIMBAD lists the spectral type as B 1.5 V. This star 
has a l so been classified as B2 V ( Stebbins, H uff'er & WhitfordI 
ll940l:ICrawfordll958l) . IWarren & Hesser,(,1977,) assign it a spec- 
tral class between B 1 .5 V to B2 V. and note that it may be vari- 
able in V, which is confirmed by examining t he me asurements 
in the the SIMBAD data base. IStephensonI (Il986i) identified 
HD 37903 as a weak Ha emissi on line star in an objecti ve prism 
survey of the northern sky. Yet ' Hernandez et alj (l2005h did not 
classify it as a Herbig Ae/Be star in their study of early type stars 
in nearby OB associations. Since HD 37903 illuminates a reflec- 
tion nebula, is an early B-star, is photometrically variable and as 
we show below, has an infrared excess, it appears to have all the 
characteristics of a Herbig Ae/Be star (iThe. de Winter & Perez! 
1994). Although it is a very bright star (V - 7.8™) and extremely 
well studied, almost all spectral classifications of the star are 
based on photometry. There are, however, several long integra- 
tion UV spectra of excellent quality available in the archives 
of the International Ultraviolet Explorer (lUE), which enabled 
us to perform an improved spectral classification. We therefore 
looked at the short wavelength high dispersion (SWP) spectra in 
the lUE ar chive and used the spectral classification of B stars de- 
veloped bv lRountree & Sonn eborn (1 99 1 ) for ultraviolet spectra 
to derive a more accurate spectral classification. This compari- 
son shows that the star is definitely later than B1.5 and earlier 
than B2.5, i.e. a B2 V star is a very good match. We therefore 
conclude that HD 37903 is a Herbig Ae/Be star of spectral class 
B2 Ve. 

To verify that the star has an infrared excess we compared the 
observed infrared fluxes to the Kurucz stellar atmosphere model 
corresponding to solar metalicity, log[g]=3.0 and an effective 
temperature of 21000 K (Fig. ^middle). We find that while 
the optical and NIR fluxes are satisfactorily reproduced by the 
Kurucz model, the IRAC fluxes from 4.5 fim show clear infrared 
excess compared to the model and a very large excess at 24 jum 
possibly indicating the presence of a protostellar disk. Although 
the star is clearly seen at 24 fj.m, the measured flux density is 
very uncertain due to the intense PDR ridge just sout h of it. 
Since HD 37903 illuminates a faint H ii region ( Wvrowski et al.l 
2000), some of the 24 //m emission also originate in the H n re- 
gion rather than from a disk. Although HD 37903 has an IR ex- 
cess, we see no trace of it in the sub-millimeter continuum. This 
is not surprising, since disks around early B-stars are likely to be 
short lived. We do not know of any early B-star, that illuminates 
a large reflection nebula a nd still would be associat ed with a cold 
protostellar disk (see e.g. lFuente et al.lll998ll200l[) . 

V615 Ori or L1630MIR-52 is a heavily reddened T Tauri 
star w ith spectral type of G7 III-IV (Malin, Ogura & Wals^ 
119871) . It is one of the few PMS stars in the field that is a gen- 
uine IRAS source (IRAS 05388-0224) and not a false detec- 
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Fig. 11. Gray-scale 850 and 450 /im continuum maps overlaid with contours of the field surrounding NGC 2023 MM 5. The 
peak intensities at 850 and 450 yum are 0.88 and 5.38 Jy/beam respectively. The contours for both the 850 and 450 fim maps 
are at 20 to 100% of the respective peak intensities in steps of 10%. The beamsizes at 450 and 85 0jum are shown at the right 
bottom corner of each panel. "+" and "x" signs mark the sub-mm sources detected in this region by iJohnstone et al.l (l2006l) and 
iNutter & Ward-ThompsonI ( |2007|) . respectively. 



tiorO, with flux densities of 1.45 Jy at 12 fim and 1.67 Jy at 
25 fim. It is listed in the 2MASS point source catalogue and 
we detect it in all IRAC ban ds and in the MIPS 24 fim image. 
iMah n. Ogura & WalshI (1198 7) quote an observed B-V color of 2 
(sugg esting Eb-v - 1.2). In the General Catalogue of Variable 
stars (lKukarninetalJll97lb the listed V magnitudes vary be- 
tween 13.8 and 15.9. These va lues agree well with more re- 
cent data of iDroege etall (l2007l) . who give V =13.11±0.26 and 
/ =10.99+0.16. The star is also an X-ray source, detected by 
the XMM-Newton telescope (XMM-Newton 2nd Incremental 
Source Catalogue). Using the optical and infrared flux densi- 
ties and the online SED fitter we find 7 models for which x^- 
^besi ^ ^-^^ However 6 out of the 7 models have stellar temper- 
atures > 6000 K, which is inconsistent with the known spectral 
type of V615 Ori. We thus present here the only model which 
corresponds to a stellar temperature agreeing with previous ob- 
servations. The parameters for this model are: Distance - 417 pc, 
M. = 4.0 Mo, = 5546 K, = 62.7 Lg, Menv = 0.056 Mg 
and Mdisk = 0.006 Mq. If we use the power-law dependence 
between disk mass and 850 ^i m flux density derived from 153 
T Tauri stars in Taurus-Auriga dAndrews & Williamsll2005l) . we 
find that the modeled disk and envelope mass would result in an 
850 fj.m flux density of ~ 4 mJy, which is below our detection 
limit. It is therefore not surprising that we did not detect the star 
in our SCUBA observations nor would it hav e been detected in 
the deeper 850 //m image by Jo hnstone et al.l ([2006). 

iMaUn. Ogura & WalshI (Il987l) found two HH objects 
(NGC 2023-HH2, and HH 3) souflieast of V615 Ori. The proper 
motions of HH2, however, suggest that the exciting source 
is south or southeast of HH2 rather than in the direction of 
V615 Ori. There is a very red MIR source, L1630MIR-51 



^ IRAS 05384-0229 and IRAS 05386-0229 in the Horsehead neb- 
ula; IRAS 05385-0222, IRAS 05384-0219, and IRAS 05385-0212 in 
the IC434 ionization front; IRAS 05392-0214, IRAS 05390-0214, and 
IRAS 05391-0217 (near- HD 37903) in NGC2023 are all spurious de- 
tections, since they are not seen in any IRAC bands nor in the MIPS 24 
image 



(Table |2]l, south-east of HH 2, which could be the exciting star 
for HH2, and possibly HH 3. L1630MIR-51 has very strong in- 
frared excess and colors of a Class II/I object, which indicate 
that it is very young and deeply embedded. Examination of the 
large CO 7 = 2^1 map to which we have access, indicates that 
V615 Ori may drive a weak, extended CO outflow. A position- 
velocity plot through the star and HH 3 shows faint blue-shifted 
emission to the south-east, i.e. in the direction of HH 3, and red- 
shifted emission to the northeast, with a reversal of the line wings 
approximately at the position of V615 Ori. Nothing stands out 
in the CO map at the position of HH2 or L1630MIR-51. We 
therefore believe that HH3 is excited by V615 Ori, although 
we caution that the identification is tentative and has to be con- 
firmed with deeper observations with better spatial resolution. 
HH2 may be part of the V615 Ori outflow, although we can- 
not exclude that this HH complex is excited by L1630MIR-56, 
which is ~ 30"-40" south east of HH2 A,B. 

Sellgren's star C (L1630 MIR-62, HBCSOtB) was als o 
found to be a T Tauri star by iMalin. Ogura & WalshI ( Il987h . 
but with much stronger Ha emission than V615 Ori. However, 
since the star is heavily obscured, they were able to only con- 
clude that the star has a spectral type of late G or early K and 
a luminosity class V , based on weak Na I and Mg I absorption. 
IWeaver & Babcoc^ (|2004|) also found it on a very deep CTIO 
objective prism plate from 1976. In their catalogue it is number 
54, and they characterize the star as having medium strength Ha 
emission with weak continuum. The star is rather bright in IRAC 
and MIPS images, L1630MIR-62, (Tabled, and has MIR colors 
of a Class II object. 

NGC 2023-VLA 3 or L1630MIR-46 was the stronge st radio 
source (~ 50 mJy) in the 3.6 cm VLA survey by Reipurth et al.l 
(120041) . It has no optical or sub-millimeter counterpart, but is 
detected in all IRAC bands and also in the 24 yum MIPS im- 
age (Table |2]l with colors of a Class II or borderline Class I 



* Both Sellgren's star C ( ISellgreiJl983h and HBC 500 have incorrect 
coordinate entries in the SIMBAD data base; neither corresponds to a 
visible star or IR source. 



16 



Mookerjea et al.: YSOs near NGC 2023 



source. The flux density is too high to be free-free emission from 
an ionized jet. It is almost certainly gyrosynchrotron emission, 
which is most often found in weak-line (naked) T Tauri stars, but 
sometimes in young stellar objects with strong magnetic fields 
jAndre et al.|[l992 ). Since the star is an X-ray source (XMM- 
Newton 2nd Incremental Source Catalogue), it is likely to be 

have magnetic activity. 

Additional sources: IWeaver & Babcockl (|2004|) identified 
seven additional Ha emission line stars within our field of view, 
WB 53 - 59. WB 53 is most likely the Herbig-Haro object HH 5 
all the others have MIR counterparts and all but two, WB57 
and WB 59, have clear IR-excess. VLA2 has no near- or mid- 
IR counterpart and is most likely an extragalactic source. 

10. Discussion 

10.1. The YSO population 

Extracting photometry of "faint" stars in a bright knotty reflec- 
tion nebula like NGC 2023 is rather involved. The Spitzer pho- 
tometry package APEX recovered far more sources than what 
we present in Table |2] However, many of the "detections" had 
blue [5.8]-[8.0] colors and very red [4.5]-[5.8] colors, i.e. they 
appeared non-stellar. Visual inspection of each source showed 
that they all appeared nebulous and fuzzy. It is quite likely that 
some of them may contain real stars. The PSF photometry, how- 
ever, was definitely dominated by nebular emission and/or ab- 
sorption present in the broad IRAC filters, and we therefore 
rejected them. The final list may still contain some spurious 
sources, but if there are, they are very few. 

With a total of 73 YSOs identified in the NGC 2023/L 1630 
the sample is sufficiently large to visualize the distribution of the 
YSOs (Fig. |2]i. We find that the majority of them, ~ 23 stars, 
are located in the southern part of NGC 2023. There is another 
7-9 young stars in the sub-millimeter ridge west of NGC 2023, 
in which the Class protostar, MM 1, is embedded. If we con- 
sider these to belong to the NGC 2023 cluster, we therefore have 
a total of 32 YSOs. We also find a sti'ing of YSOs (afl Class II) 
along the whole submillimeter ridge, which runs parallel to the 
bright FDR rim IC 434. A sUghtly larger fraction of the YSOs 
lie on the western side of this dust ridge, i.e. the side facing the 
large IC 434 H ii region. The rest of the YSOs lie in the bright 
rim created by IC 434 or slightly inside IC 434. Here we observe 
two distinct groups. One group is associated with the Horsehead 
nebula, where we see 7 sources just outside the head region and 
another two down the neck, where the Horsehead nebula merges 
into the L 1630 molecular cloud. Another small group is seen in 
the northwestern part of our field. The latter group is centered 
around an eroded elephant trunk structure, which is now much 
smaller than the Horsehead. A few PMS stars (including VLA 3) 
are scattered all over the L 1630 dark cloud, suggesting that oc- 
casional star formation can occur in any dense cold cloud core 
without any external triggering. Of the 5 bona fide Class I/O stars 
from Table|2]and the 10 very young MlPS-only sources, a total 
of 9 are concentrated in the southern part of NGC 2023 and near 
or in the MM 1 cloud core. The remaining 6 are distributed in a 
scattered fashion close to the Horsehead nebula and the T Tauri 
star V615 Ori. The strong concentration of Class I/O sources near 
NGC 2023 suggests that this is where the youngest stars are and 
where star formation may still be on-going. 

Our YSO sample includes all the Ha emission line objects 
from the Weaver & Babcock (2004) survey. However the two 
emission line objects, WB 57 and WB 59, which are inside the 
NGC 2023 reflection nebula and are to faint to be detected at 



8.0 yum, but are readily detected at 3.6 and 4.5 jum. Since both lie 
in areas of intense PAH emission, it is possible that the neutral or 
slightly blue [3.6] - [4.5] color is due to intense PAH emission 
in the 3.6 yum band. As we have no way of confirming this at 
the present, we have counted them as members of our candidate 
YSO population. 

Our study re-enforces the view that a large H ii region like 
IC 434 expanding into dense surrounding cloud complexes trig- 
gers the formation of new stars. The densest regions in the sur- 
rounding molecular clouds, like the Horsehead nebula, can even 
form small clusters. The Horsehead may even form additional 
stars, because we found one Class I/O source near the base of 
the Horsehead. There is also a bright, extended pre-stellar core 
in the neck of the Horsehead, which is seen as an infrared dark 
cloud in IRAC and MIPS images. The pre-shock from the Hii 
region can compress molecular gas even further into the cloud 
and trigger another wave of star formation. This is demonstrated 
by all the young stars we find in submillimeter ridge, which runs 
parallel to the IC 434 boundary, but almost 0.8 pc from the edge 
ofIC434. 

A comparison of our results with the recent study to de- 
rive an infrared census of YSOs within a 7'x7' region around 
the Horsehead nebula dBowler et alj|2009l) . shows a very good 
agreement with the results we get from IRAC photometry. 
We h ave identified almost all the sources with an IR excess, 
whichl Bowler et al] ( |2009l) present in their Table 1. The sources 
which we do not detect are B33-13, B33-20,B33-27, B33-29, 
and B33-30 all of which were selected on the basis of their 
NIR colors, a criterion which does not match with our MIR- 
bas ed approach. Th e classification of these sources derived by 
Bowler et a 120091) also match with the classes derived in this 
paper. Although iBowler et alj (l2009h report a significant 8 jum 
excess for B33-31 and classify it as a bona-fide YSO, we do not 
detect the source at all at 8.0 jum. We do retain it as a candidate 
YSO. Since we have also analyzed MIPS data for the region, 
we do identify a lar ger number of candidate YSOs than seen by 
iBowleretalJ ( 120091) . 

The present Spitzer-hased study has enhanced the number of 
YSOs in the young star cluster associated wi th NGC 2023 , which 
has HD 37903 as its most massive member. iDepov et alj (Il990h 
found 16 young stars in this cluster based on deep near-IR imag- 
ing. We find that two of them are heavily reddened foreground 
stars, and the true number is therefore 14. In this study we find 
~ 32 or possibly as few as 22 YSOs in the star cluster associated 
with NGC 2023. Even with our larger estimate of the number 
of YSOs in the NG C 2023 cluster, we are not even close to the 
number of stars, 50. lDepov et al.l(ll990h had estimated to be sur- 
rounding HD 37903. Their estimate was based on the assump- 
tion that even a small cluster foll ows a Salpeter or Mille r & Scalo 
IMF, intersteUar mass function dMiller & Scalolll979l) . Several 
investigations, especially from groups that study clusters in star 
forming galaxies, show that there can be substantial stochastic 
fluct uations when one deals with small number sta tistics (see 
e.g.. lCervino & VaUs-Gabaudll2003l: ICarigi & Hernandez.20081) . 
However, we do know that we have missed PMS stars at the faint 
end, i.e. very low mass stars. If we add stars, which were not 
seen by 2MASS, but which are very red (> 0.3) based on their 
[3.6]-[4.5] colors, we find another 8 stars south of the reflection 
nebula, one to the north west of HD 37903 inside the reflection 
nebula, and another 7 in the sub-millimeter ridge. Most of these 
are likely to be low mass PMS stars. The true, but still conserva- 
tive number for the NGC 2023 cluster is therefore ~ 32 stars. If 
we account for stochastic fluctuations and also take into consid- 
eration that this a very young cluster, where star formation is still 
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going on, there is nothing anomalous with the mass distribution 
of these young stars. Some of the low-mass stars may still be in 
the process of forming or be deeply embedded in the molecular 
cloud core. They may therefore not even be detectable by deep 
mid-lR observations, like the IRAC observations that we have 
employed in this study. 

10.2. The sub-millimeter and protostellar candidates 

SCUBA maps at 450 and 850 fum reliably detect only 5 sub- 
millimeter sources, MM 1-5. Of these sub-millimeter cores, 
MM 1 has been shown to be a Class object by ISandell et alJ 
(11999.) . Based on isothermal grey-body modelling we find that 
MM 1 has a rather low value of /?. Since MM 1 is known to drive 
an energetic outflow, it must have an actively accreting disk. It 
therefore seems that the inferred low p value is unlikely to be 
due to the formation of large grains. We argue that the low /3 
arises due to the inner disk being partially optically thick. We 
have imaged the protostellar disk associated with MM 1 using 
interferometric millimeter observations of the NH2D emission. 
However based on the observed velocity structure we conclude 
that the NH2D emission appears to trace rotation in the surround- 
ing envelope rather than the protostellar disk. The MM 1/MM 2 
region coincides with an IR dark cloud and the presence of two 
faint MIPS sources suggests that stars have already formed in the 
cloud core. MM 2 and the surrounding cloud core is extremely 
dense and cold and is not evident in single dish C'^O and C'^O, 
while NH2D shows a very similar morphology to that of the dust 
continuum, especially for the cold, dense, dark cloud core sur- 
rounding MM 1 and MM 2. Analysis of column densities derived 
from C'**0 and C'^O J = 2 ^ 1 compared to column densities 
estimated from the dust emission suggest that CO is strongly de- 
pleted in the cold core, presumably because a large fraction of 
the CO is frozen onto dust grains. NH2D, on the other hand, ap- 
pears to be strongly enhanced in the cold cloud core. The mass 
of the gas of MM 2 and the cloud core surrounding MM 1 and 
MM 2, ~ 25 Mo, is certainly sufficient to form additional low- 
mass stars in this cloud core. 

The sub-millimeter sources MM 3 and MM 4 lie next to the 
NGC 2023 reflection nebula, which is where we see the high- 
est concentration of young stars. Here the surrounding molec- 
ular cloud is much warmer and we detect no CO depletion. 
MM 3 is a compact sub-millimeter source, which coincides with 
a bright mid-infrared source, L1630MIR-63, which has a very 
strong infrared excess and which we classify as an extreme Class 
I source based on color-color diagrams. In the millimeter and 
sub-millimeter is looks almost identical to MMl, although it 
is not clear whether it drives an outflow and it does not have 
any observable free-free emission. The dust emissivity index, /3 
= 0.8, deduced for MM 3 is almost identical to MM 1 and al- 
though the disk/envelope mass is only about a third of the disk 
around MM 1, both stars have very massive disk. Yet MM 3 is 
detected even in the near-lR, while MM 1 is not even seen at 24 
fim, suggesting that MM 3 is a more evolved star than MM 1 . 
It is, however, possible that the main reason why we do not de- 
tect MM 1 in the mid-IR is because we see the disk nearly edge 
on. If the disk surrounding MM 3 is more edge on, and there- 
fore not blocking the light from the central protostar, both could 
be in a very similar stage of their evolution. MM 4 is the most 
extended dust core we identified in our sample. It differs from 
typical pre-stellar cores, because it is relatively warm and it also 
has a slightly anomalous y6-index. MM 5 is a cold prestellar core 
with both temperature and /3 characteristic of most pre-stellar 
cores. 



As discussed briefly in Sees. |5] and [T] we find 
substantially different masses for almost all o f the sub- 
millim eter c ores that ove rlap with the sample of John stone et al.l 
(120061) and iNutter & W ard-Thompson (2007), except MM 5, 
(SMM J05414 9-02212 (Johnstone et al. 2006), and OrionBS- 
541299-221 14 dNutter & Ward-Thompson 2007)). For this sub- 
millimeter core all three studies find about the same mass, ~ 
3-4- Mq. Here we review why there are such large difference in 
masses estimated for the cores MM 1- 4. 



iJohnstone et a D i2006) and iNutter & Ward-ThompsonI 
( 12007 ) estimate very similar masses for the MM 1/MM 2 re gion, 
~ 15 Mo (SMM J054141-02180 iJ^ohnstone et aTl I2OO6I) and 
~ 13 Mq (OrionBS-541257-21809 iNutter & Ward-ThompsonI 
I 2OO7I) , while we show that this "core" breaks up into two 
sub-millimeter cores embedded in a very dense, cold cloud. 
MM 1 coincides with a Class protostar, while MM 2 appears 
to be a pre-stellar core (Section |5j. They diff'er in temperature 
and have very diff'erent SEDs. MM 2 has a /? ~ 2, while MM 1 
has a very flat SED with /? - 0.7. If we add the masses for MM 1 
and MM 2 we get ~ 8 Mo, only about half of the mass found 
by treating MM 1 and MM 2 as a single core. For MM 3, which 
we also identify as a Class O/I protostar, with an SED similar 
to M M 1 , the differences in mass are large. We find 0.3 M©, 
while Hohnsto ne et al. (2006) determ ine a mass for of ~ 16 Mp 
core SMM J054 162-02 172), and INutter & Ward-ThompsonI 
Soot!) find 8 Mo (core OrionBS-541375-21733). Here the large 
diff'erence is in part due to MM 3 having a very flat SED, but 
also because both groups deduce a much larger size for MM 3 
than we do. We find it point-like and unresolved (Table [TJ. 
For MM 4 Johnstone et al. (2006) (SMM J054167-02171) and 
Nutter & Ward-Thompsom (,2007.) (OrionBS-541404-21702) 
obtained masses of 8 and II.OMq, respectively, which is about 
10 times higher than what we find. In this case the discrepancy 
primarily arises due to the differences in the derived source 
size and assuming normal dust properties, while we find that 
the dust emissivity diff'ers from that of a typical pre-stellar core 
(SectionO. 



We attribute the primary source of the different mass esti- 
mates to the manner in which a core is defined. Automatic clump 
finding algorithms like clfind appear to find large cores, if the 
core is embedded in extended emission. Furthermore the previ- 
ous studies which are both at a single wavelength also identify 
"clumps", which appear spurious, like some of the cores in the 
molecular cloud south of NGC 2023. A protostellar core should 
have a correspondence in some dense gas tracers, like HCN, 
HCO^, CS, etc. and stand out with a single velocity. It can have 
a gradient due to collapse, rotation or expansion. All the cores 
do not necessarily have the same temperature or dust properties. 
Cores with embedded protostars have rather different dust prop- 
erties (j6-index) and temperature structure. Resolution of the ap- 
parent disagreement between different sets of observations will 
be possible only by using far-infrared observations having simi- 
lar spatial resolution, so that it is possible to determine the dust 
emissivity and temperatures of these cores reliably. This will 
soon be possible with missions like Herschel and SOFIA, both of 
which have far-infrared imagers providing spatial resolution at a 
100 fim of 8 - 10", very similar to what can be achieved on large 
ground based single dish telescopes. Only then can we reliably 
assess how closely the pre-stellar core mass function resembles 
the interstellar mass function. 
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11. Summary 

We have investigated the distribution of gas, dust and young stars 
in the vicinity of NGC 2023 and the L 1630 molecular cloud 
using a combination of the best-available datasets at multiple 
wavelengths. 

We reliably identify 5 sub-millimeter cores MM 1-5 at 450 
and 850 fim using SCUBA observations. Of these cores MM 1 
and MM 3 are bona fide Class 1/0 sources. Single tempera- 
ture greybody models have been fitted to the SEDs of the sub- 
millimeter cores to derive dust temperatures and emissivity in- 
dices (J3). We find that adjacent cores MM 1 and MM 2 show 
vastly different /3 values of 0.7 and 2.1 respectively. Large dif- 
ferences in jS values are obtained for the sources MM 3 (0.8) and 
MM 4 (1.5) as well. The lower (< 1) values of /3 compared to 
those found in PMS objects could be explained by the presence 
of partially optically thick disks and/or grain growth in disks. In 
either case it is a clear signature of the presence of protostellar 
disks. We find that the values of /3 correlate with the dust tem- 
peratures and hence the evolutionary stage of the source. The 
derived significant variation of /3 values highlights a potential 
problem in the commonly made assumption of uniform values 
of jS when calculating masses of cores distributed over a region. 

Using NIR (2MASS) and MIR (Spitzer IRAC and MIPS) 
photometry we have generated color-color diagrams to estimate 
the ages of the detected sources. We identify a total of 73 PMS 
objects in the entire 22' x 26' region, out of which a possible 
22-32 are associated with the very young cluster associated with 
NGC 2023 of which HD 37903 is the most massive member We 
categorize 5 YSOs as Class 0/1 objects and 10 sources as Class 
I/ll objects. For the 10 sources which are detected only in the 
MIPS 24 jim band the protostellar nature needs to be confirmed. 
Except for MM 3 none of the young MIR sources has been de- 
tected in the sub-millimeter. This, together with the fact that the 
available low (compared to the present datasets) resolution FIR 
maps look significantly different from the data presented here, 
emphasizes the need for follow-up FIR observations with com- 
parable resolution, in order to conclusively determine the evo- 
lutionary stages of these objects. The MIR datasets also reveal 
several sources like HD37903, V615 Ori, Sellgren's star C etc., 
which are bright in the optical. The UV datasets have been used 
to refine the spectral classification of HD37903 and the MIR 
datasets further suggest that it as a Herbig Be star with signif- 
icant infrared excess. For V615 Ori, which most likely drives a 
weak extended CO outflow, the infrared excess can be explained 
by a low-mass (0.06 Mq) circumstellar disk and envelope, which 
is too faint to be detected in the sub-millimeter. 

We have used high resolution NH2D observations to probe 
the protostellar disk associated with the Class source MM 1 . 
Although we detect enhanced NH2D emission, the velocities 
suggest that the emission arises primarily from the surrounding 
ambient cloud which is also rotating. Based on C'^O observa- 
tions we confirm CO depletion in the MM 1/MM2 region, con- 
sistent with the enhanced NH2D emission. In contrast MM 3, 
also identified to be a Class O/I star and estimated to have 
a dust temperature, somewhat higher than the temperature of 
MM 1/MM 2 shows no signs of CO depletion. 

Based on the derived distribution of YSOs in 
NGC 2023/L 1630 region we conclude that the expansion 
of the H II region IC 434 has triggered the formation of stars in 
the compressed gas inside the L 1630 molecular cloud, as well 
as in some of the dense "elephant trunk" structures. We have 
been able to derive reliable age estimates of most of the YSOs in 
the mapped region based on NIR and MIR photometry. However 



far-infrared observations at comparable angular resolutions are 
absolutely essential to determine the evolutionary stages of the 
extremely young objects which are clearly detected only at 
24 fim and also to confirm the protostellar nature of three of the 
sub-miflimeter clumps (MM 2, MM 4 and MM 5) for which no 
other signatures are so far available. 
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Table 2. Coordinates and flux densities of MIR sources detected in the vicinity of NGC 2023 



Source 


<j;2ooo 


^2000 


/ 


H 


Ks 


F3.6 jmi 


F4.5 pan 


F5.8 prni 


Fs.o Jim 




Other Name 


Class 


L1630 












mJy 


mJy 


mJy 


mJy 


mJy 






MIR-1 


5:40:39.31 


-2:26:45.8 


11.58 ±0.03 


10.70 ± 0.02 


10.31 ± 0.02 


28.71±0.10 


26.48 ± 0.07 


31.0 ±0.6 


46.4 ± 0.4 


70.2 ± 0.2 




n 


MIR-2 


5:40:40.01 


-2:12:46.7 


12.20 ± 0.03 


11.45 ± 0.02 


11.23 ±0.02 


10.41±0.06 


6.14 ±0.03 


5.3 ± 0.6 


2.3 ± 0.2 






s 


MIR-3 


5:40:41.39 


-2:29:07.8 


13.07 ± 0.03 


12.44 ± 0.03 


12.15 ± 0.03 


5.45±0.02 


3.71 ± 0.02 


3.8 ± 0.5 


1.3 ± 0.2 






s 


MIR-4 


5:40:41.42 


-2:19:11.2 


10.29 ± 0.03 


9.61 ± 0.03 


9.40 ± 0.02 


51.87±0.14 


32.27 ± 0.09 


23.2 ± 0.7 


13.7 ± 0.3 






s 


MIR-5 


5:40:44.56 


-2:30:36.5 


11.95 ±0.03 


11.52 ± 0.02 


11.36 ±0.02 


8.19±0.03 


5.15 ±0.02 


4.2 ± 0.4 


2.2 ± 0.2 






s 


MIR-6 


5:40:44.63 


-2:22:43.9 


11.82 ±0.03 


11.14 ± 0.02 


10.91 ± 0.03 


13.17±0.05 


8.73 ± 0.03 


7.8 ± 0.4 


3.4 ± 0.2 






s 


MIR-7 


5:40:44.75 


-2:11:49.8 


9.85 ± 0.03 


9.20 ± 0.02 


8.96 ± 0.03 


117.20±0.56 


102.20 ± 0.51 


104.5 ±1.0 


112.4 ± 0.7 


74.2 ± 0.2 




n 


MIR-8 


5:40:44.75 


-2:13:53.7 


11.49 ±0.03 


10.64 ± 0.03 


10.39 ± 0.02 


21.82±0.07 


12.76 ± 0.04 


11.0 ±0.4 


5.5 ± 0.2 






s 


MIR-9 


5:40:44.85 


-2:12:15.6 


9.76 ± 0.03 


9.00 ± 0.02 


8.80 ± 0.03 


88.02±0.28 


56.26 ± 0.15 


42.8 ± 0.7 


24.1 ± 0.4 


3.2 ± 0.2 




s 


MIR-10 


5:40:46.07 


-2:28:03.3 


12.08 ± 0.03 


11.42 ± 0.02 


11.24 ± 0.03 


9.46±0.03 


5.74 ± 0.02 


4.4 ± 0.4 


3.1 ± 0.2 






n/m 


MIR-1 1 


5:40:46.78 


-2:10:49.9 


11.93 ±0.03 


11.02 ±0.03 


10.51 ± 0.03 


30.03 ±0.11 


27.24 ± 0.07 


23.9 ± 0.6 


26.3 ± 0.3 


72.2 ± 0.2 




n 


MIR- 12 


5:40:47.07 


-2:24:48.6 


11.51 ±0.03 


11.17 ± 0.03 


11.10 ±0.03 


10.53±0.03 


6.97 ± 0.03 


5.4 ± 0.4 


2.1 ±0.2 






s 


MIR-13 


5:40:47.53 


-2:12:10.9 


12.36 ± 0.03 


11.58 ±0.02 


11.36 ±0.03 


8.40±0.04 


5.34 ± 0.02 


4.5 ± 0.5 


2.1 ± 0.2 






s 


MIR-14 


5:40:49.61 


-2:28:56.0 


9.70 ± 0.03 


8.98 ± 0.02 


8.70 ± 0.02 


92.42±0.29 


58.86 ± 0.13 


42.7 ± 0.7 


24.0 ± 0.3 


3.2 ± 0.2 




s 


MIR-15 


5:40:50.48 


-2:12:53.9 


9.88 ± 0.03 


8.99 ± 0.02 


8.74 ± 0.02 


97.28±0.52 


65.38 ± 0.41 


46.5 ± 0.7 


22.7 ± 0.4 






s 


MIR-16 


5:40:50.94 


-2:09:45.2 


11.44 ±0.03 


10.76 ± 0.02 


10.53 ± 0.02 


20.54±0.07 


12.94 ± 0.04 


10.3 ± 0.5 


8.4 ± 0.3 






n/m 


MIR- 17 


5:40:51.10 


-2:20:59.8 


11.30 ±0.03 


11.02 ±0.02 


10.95 ± 0.03 


10.93±0.04 


7.49 ± 0.02 


6.2 ± 0.5 


2.9 ± 0.2 






s 


MIR- 18 


5:40:51.73 


-2:26:48.9 


12.25 ± 0.04 


10.73 ± 0.03 


9.71 ± 0.03 


142.00±0.65 


166.40 ± 0.58 


217.4 ± 1.4 


322.3 ±1.1 


754.5 ± 0.5 


B33-1 


n 


MIR- 19 


5:40:52.41 


-2:27:12.6 


15.09 ± 0.05 


13.61 ± 0.04 


12.58 ± 0.03 


7.33±0.04 


9.12 ±0.03 


12.7 ± 0.5 


22.2 ± 0.3 


105.0 ± 0.2 


B33-28 


W 


MIR-20 


5:40:53.11 


-2:15:48.4 


11.02 ±0.03 


10.51 ± 0.02 


10.29 ± 0.02 


23.95±0.07 


15.91 ± 0.05 


11.9 ± 0.4 


6.1 ± 0.2 






s 


MIR-21 


5:40:53.75 


-2:09:38.4 


12.27 ± 0.03 


11.29 ±0.02 


10.86 ± 0.02 


15.85±0.05 


11.40 ±0.03 


7.9 ± 0.4 


3.1 ± 0.2 






s 


MIR-22 


5:40:54.86 


-2:24:06.6 


12.04 ± 0.03 


11.27 ± 0.02 


11.06 ± 0.02 


11. 00 ±0.04 


7.28 ± 0.03 


6.6 ± 0.5 


3.1 ±0.2 






s 


MIR-23 


5:40:56.09 


-2:24:02.1 


9.47 ± 0.03 


8.79 ± 0.05 


8.57 ± 0.02 


101.30±0.52 


71.04 ± 0.35 


48.1 ± 0.6 


27.3 ± 0.3 


1.30±0.1 




s 


MIR-24 


5:40:56.68 


-2:11:15.7 


12.67 ± 0.03 


10.88 ± 0.02 


10.07 ± 0.02 


43.18±0.14 


29.22 ± 0.09 


22.7 ± 0.5 


15.3 ± 0.4 






n/m 


MIR-25 


5:40:56.73 


-2:08:57.1 


8.91 ±0.02 


8.74 ± 0.04 


8.64 ± 0.02 


97.09 ±0.51 


64.84 ± 0.36 


42.4 ± 0.7 


21.2 ± 0.4 




HD 2908 15 


s 


MIR-26 


5:40:56.89 


-2:24:49.6 


13.22 ± 0.03 


12.51 ± 0.02 


12.27 ± 0.02 


3.49±0.02 


2.27 ±0.11 


1.9 ±0.5 


1.0 ± 0.2 






n/m 


MIR-27 


5:40:58.14 


-2:31:32.2 


12.38 ± 0.03 


11.64 ±0.03 


11.46 ±0.02 


8.41±0.03 


5.20 ± 0.02 


4.3 ± 0.4 


2.1 ± 0.2 






s 


MIR-28 


5:40:58.68 


-2:25:26.7 


12.25 ± 0.03 


10.77 ± 0.02 


10.06 ± 0.02 


39.42±0.14 


30.27 ± 0.09 


21.2 ± 0.6 


13.5 ± 0.3 




B33-10 


n/m 


MIR-29 


5:41:01.57 


-2:21:13.9 


10.47 ± 0.03 


8.81 ± 0.03 


8.05 ± 0.03 


218.80±1.00 


143.40 ± 0.64 


113.7 ±0.9 


63.6 ± 0.6 






s 


MIR-30 


5:41:01.92 


-2:07:03.4 


13.96 ± 0.03 


12.21 ± 0.02 


11.34 ±0.02 


13.30±0.06 


10.18 ± 0.04 


7.4 ± 0.5 


3.5 ± 0.2 






n/m 


MIR-3 1 


5:41:02.67 


-2:18:17.8 


6.95 ± 0.03 


6.87 ± 0.05 


6.76 ± 0.02 


517.60±1.86 


337.60 ±1.18 


228.3 ± 1.2 


131.0 ±0.7 


55.4 ± 0.2 


HD 37805 


s 


MIR-32 


5:41:03.57 


-2:25:12.5 


11.81 ±0.03 


10.86 ± 0.03 


10.38 ± 0.02 


23.72±0.08 


16.13 ± 0.05 


11.8 ±0.5 


7.7 ± 0.3 






s 


MIR-33 


5:41:03.78 


-2:06:26.3 


12.85 ± 0.03 


12.25 ± 0.02 


11.94 ± 0.02 


7.30±0.03 


5.90 ± 0.03 


5.7 ± 0.4 


4.8 ± 0.3 


15.0 ± 0.2 




n 


MIR-34 


5:41:05.44 


-2:31:05.9 


12.74 ± 0.03 


11.75 ± 0.03 


11.40 ±0.02 


9.72±0.03 


5.90 ± 0.02 


6.0 ± 0.4 


6.9 ± 0.3 






n/m 


MIR-35 


5:41:05.44 


-2:24:59.0 












1.3± 0.4 


3.2 ± 0.3 






n? 


MIR-36 


5:41:09.41 


-2:14:34.0 


11.79 ±0.03 


10.96 ± 0.02 


10.60 ± 0.02 


19.74±0.06 


13.60 ± 0.05 


11.1 ±0.5 


5.5 ± 0.3 






s 


MIR-37 


5:41:10.31 


-2:29:04.0 


13.83 ± 0.03 


13.27±0.03 


12.91±0.03 


2.39±0.02 


1.83±0.01 


1.27±0.03 


1.22±0.2 




B33-21 


n/m 


MIR-38 


5:41:12.23 


-2:27:34.5 


13.70 ± 0.03 


12.12 ±0.03 


11.39 ± 0.03 


12.08±0.05 


7.77 ± 0.03 


6.8 ± 0.4 


4.6 ± 0.2 




B33-22 


n/m 


MIR-39 


5:41:12.40 


-2:19:21.1 








0.35±0.01 


0.35 ± 0.01 


1.4 ± 0.4 


2.6 ± 0.2 






W 


MIR-40 


5:41:13.06 


-2:29:12.2 


11.78 ±0.03 


11.37 ± 0.02 


11.21 ±0.02 


9.22±0.03 


6.42 ± 0.02 


4.7 ± 0.4 


2.2 ± 0.2 




B33-23 


s 


MIR-4 1 


5:41:14.06 


-2:25:05.1 


12.76 ± 0.03 


12.11 ± 0.03 


11.91 ± 0.02 


5.78±0.24 


3.57 ± 0.02 


3.4 ±0.1 


4.1 ± 0.2 




B33-25 


n 


MIR-42 


5:41:14.45 


-2:29:35.4 


10.81 ± 0.03 


10.51 ± 0.03 


10.43 ± 0.02 


17.24±0.06 


12.20 ± 0.04 


8.7 ± 0.5 


3.9 ± 0.2 






s 


MIR-43 


5:41:17.18 


-2:18:07.3 








0.35±0.00 


0.67 ± 0.01 


2.2 ± 0.4 


3.0 ± 0.2 


14.9 ± 0.2 




W 


MIR-44 


5:41:17.33 


-2:31:38.1 


9.64 ± 0.03 


9.32 ± 0.03 


9.18 ± 0.02 


60.57±0.19 


38.14 ± 0.10 


25.6 ± 0.6 


13.8 ± 0.3 






s 


MIR-45 


5:41:20.89 


-2:17:52.8 


15.43 ± 0.05 


13.35 ± 0.02 


12.16 ± 0.02 


10.61±0.04 


12.62 ± 0.04 


13.0 ±0.5 


13.8 ± 0.3 


2.0 ± 1.0 




n/m 


MIR-46 


5:41:21.70 


-2:11:08.2 








3.53±0.02 


4.84 ± 0.21 


7.8 ± 0.4 


9.1 ± 0.3 


23.4 ± 0.2 


VLA3 


n 


JVllK-4/ 


J.41.ZZ.14 


-z.io.44.U 


1 _i_ r\ f\A 
IJ.Dz ± U.UO 




1 yiQ _i_ n ni 
1Z.4(5 ± U.Uj 


cS.jU±U.Jj 


Q CO _i_ n ri'j 


Q 'J _i_ n < 
y.j ± U.J 


0.4 ± V.D 


o f\ _i_ 1 n 




11/111 


MIR-48 


5:41:22.48 


-2:17:18.6 


14.68 ± 0.04 


12.95 ± 0.03 


12.05 ± 0.02 


8.05±0.32 


6.89 ± 0.03 


5.3 ± 0.4 


2.8 ± 0.2 






n/m 


MIR-49 


5:41:22.49 


-2:09:14.7 


16.93 ± . . . 


15.74 ±0.13 


14.60 ±0.11 


0.74±0.01 


0.69 ± 0.01 


1.2 ±0.5 


3.8 ± 0.2 


13.2 ±0.1 




W 


MIR-50 


5:41:23.29 


-2:17:35.7 


17.12 ± ... 


16.53 ± ... 


14.25 ± 0.07 


4.57±0.02 


8.59 ± 0.03 


14.8 ±0.5 


19.8 ±0.4 


122.6 ± 0.2 




I 


MIR-5 1 


5:41:24.21 


-2:16:06.5 








2.61±0.01 


5.63 + 0.02 


5.3 ± 0.4 


5.0 ± 0.3 


369.8 ± 0.3 




1/11 


MIR-52 


5:41:24.49 


-2:22:36.1 


9.10 ±0.02 


8.07 ± 0.04 


7.39 ± 0.02 


922.90±2.29 


900.90 ± 2.44 


980.1 ±2.8 


897.9 ± 3.0 


1067.0 ± 0.5 


V615 Ori 


n 


MIR-53 


5:41:26.20 


-2:10:07.5 


10.46 ± 0.03 


10.22 ± 0.02 


10.13 ± 0.02 


22.49±0.08 


6.20 ± 0.03 


10.7 ± 0.5 


6.3 ± 0.2 






s 



Table 2. continued. 



Source 


12000 


52000 


J 


H 




F3.6 fim 


F4.5 /jm 


F5.8 /im 


Fg.o /jm 


F24 //m Other Name 


Class 


LI 630 












mJy 


mJy 


mJy 


mJy 


mJy 





MIR-54 
MIR-55 
MIR-56 
MIR-57 
MIR-58 
MIR-59 
MIR-60 
MIR-61 
MIR-62 
MIR-63 
MIR-64 
MIR-65 
MIR-66 
MIR-67 
MIR-68 
MIR-69 
MIR-70 
MIR-71 
MIR-72 
MIR-73 
MIR-74 
MIR-75 
MIR-76 
MIR-77 
MIR-78 
MIR-79 
MIR-80 
MIR-81 
MIR-82 
MIR-83 
MIR-84 
MIR-85 
MIR-86 
MIR-87 
MIR-88 
MIR-89 
MIR-90 
MIR-91 
MIR-92 
MIR-93 
MIR-94 
MIR-95 



5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 
5:41 



26.23 
27.98 
28.95 
30.24 
32.95 
33.90 
34.72 
36.01 
36.39 
37.18 
37.46 
37.90 
38.40 
40.25 
42.23 
42.25 
42.39 
43.10 
44.23 
44.62 
44.69 
44.79 
45.50 
45.69 
46.08 
46.68 
47.04 
48.32 
50.35 
51.40 
51.60 
52.32 
55.57 
56.32 
57.94 
59.30 
59.90 



5:42:00.01 
5:42:02.31 
5:42:02.34 
5:42:02.64 
5:42:03.77 



-2:16:22.3 
-2:19:40.2 
-2:23:19.2 
-2:07:31.8 
-2:28:02.8 
-2:19:00.3 
-2:17:24.0 
-2:08:21.1 
-2:16:46.1 
-2:17:17.2 
-2:30:21.1 
-2:25:27.2 
-2:15:32.7 
■15:59.7 
■16:24.5 
■17:35.9 
■17:21.0 
■10:26.1 
■16:16.3 
■16:06.8 
:16:56.6 
■15:55.1 
-2:24:16.0 
-2:24:36.3 
-2:31:42.8 
-2:32:05.1 
-2:16:37.8 
-2:10:40.4 
-2:31:00.8 
-2:29:51.1 
-2:07:41.3 
-2:18:32.5 
-2:23:40.3 
-2:12:45.6 
-2:11:24.2 
-2:17:02.1 
-2:21:08.7 
-2:31:48.8 
-2:23:15.6 
-2:23:28.4 
-2:07:45.8 
-2:24:50.4 



11.34 ±0.02 10.94 ±0.02 10.89 ± 0.02 



13.81 ±0.03 
15.09 ± 0.03 
18.61 ± ... 
12.28 ± 0.03 

14.76 ± 0.04 
11.73 ±0.03 

16.94 ± ... 

14.15 ±0.03 
7.37 ± 0.02 
11.31 ±0.03 
15.96 ±0.11 
12.31 ±0.02 
14.57 ± 0.03 
11.25 ±0.02 
12.86 ± 0.02 
15.19 ±0.06 

13.21 ±0.03 
13.75 ±0.08 
13.19 ±0.03 
12.61 ±0.02 

12.77 ± 0.02 
13.53 ±0.02 

9.04 ± 0.02 

10.89 ± 0.03 
13.09 ±0.02 

8.90 ± 0.02 
14.01 ± ... 

12.90 ± 0.02 
12.85 ±0.02 
12.31 ±0.02 
12.09 ± 0.03 

13.22 ±0.02 
9.99 ± 0.02 

12.22 ± 0.02 
10.67 ± 0.02 
12.70 ± 0.02 



13.18 : 
12.95 : 
15.43 : 

11.42 : 
13.38 : 

10.05 : 
13.86: 

11.46: 
7.42 : 
10.36: 
12.84: 
11.48 : 
13.20: 
10.63 : 
12.04: 

12.43 : 
12.12 : 
11.44: 
11.52 : 
11.46: 
11.64: 
12.45 : 

8.14: 
10.62 : 
12.15 : 

8.60: 
13.24: 
11.05 : 
11.66: 
11.15 : 
10.61 : 
11.67 : 

9.00 : 
11.23 : 
10.37 : 
11.88 : 



0.03 
0.03 
0.08 
0.02 
0.03 
0.02 
0.07 

0.02 
0.05 
0.02 
0.07 
0.02 
0.02 
0.02 
0.02 
0.08 
0.03 
0.05 
0.02 
0.02 
0.02 
0.02 
0.04 
0.02 
0.03 
0.02 
0.04 
0.02 
0.02 
0.02 
0.02 
0.02 
0.03 
0.02 
0.02 
0.02 



12.89 ± 0.03 
11.92 ±0.02 
13.09 ±0.03 
10.97 ± 0.02 
12.73 ± 0.02 
8.83 ± 0.02 
10.69 ± . . . 



10.18 
7.28 
9.89 
10.91 
10.99 
12.47 
10.49 
11.50 
10.38 
11.54 
9.99 
10.80 
10.87 
11.20 
11.91 
7.55 
10.58 
11.78 
8.45 
12.52 
10.16 
11.05 
10.61 
10.03 
11.02 
8.56 
10.88 
10.28 
11.48 



±0.02 
±0.02 
±0.02 
± 0.04 
± 0.02 
±0.02 
±0.02 
±0.02 
±0.03 
± 0.03 
± 0.03 
± 0.02 
±0.02 
±0.02 
±0.02 
± 0.03 
± 0.02 
±0.02 
±0.02 
± 0.04 
±0.02 
±0.02 
±0.02 
± 0.02 
±0.02 
±0.02 
±0.02 
±0.02 
±0.02 



13.15 ±0.02 11.81 ±0.02 11.28 ±0.02 



11.45±0.05 
0.27±0.00 
2.42±0.01 
9.80±0.04 
5.73±0.02 
19.50±0.08 
3.10±0.02 
505.40±1.87 
212.50±0.77 
0.22 ± 0.04 
45.88±0.14 
304.40±1.29 
23.35±0.18 
48.04±0.19 
23.17±0.09 
6.06±0.03 
19.05±0.07 
12.99±0.08 
311.00±1.00 
9.52±0.05 
137.90±0.53 
18.74±0.06 
16.69±0.05 
10.83±0.04 
7.61±0.03 
691.70±2.12 
15.52±0.06 
5.81±0.03 
131.60±0.61 
4.12±0.02 
36.53±0.11 
25.95±0.10 
18.88±0.06 
31.79±0.07 
15.61±0.06 
124.20±0.49 
14.55±0.06 
21.45±0.06 
8.62±0.03 
0.51 ±0.00 
11.01±0.04 



7.88 
1.14 
1.94 
7.44 
5.52 
17.89 
2.32 
622.50 
520.50 
0.22 
33.72 
241.20 
41.00 
56.38 
23.70 
6.20 
11.70 
11.65 
537.50 
7.96 
231.20 
13.95 
11.87 
6.87 
6.55 
721.70 
10.53 
2.69 
85.62 
2.91 
24.32 
27.56 
12.86 
24.88 
9.45 
82.14 
9.45 
13.84 
5.80 
2.88 
7.00 



±0.03 
± 0.01 
±0.01 
± 0.03 
± 0.02 
±0.06 
± 0.01 
± 1.66 
± 1.60 
±0.04 
± 0.09 
± 0.78 
±0.15 
±0.17 
± 0.07 
±0.03 
± 0.03 
±0.05 
± 1.61 
± 0.04 
± 0.89 
±0.04 
± 0.04 
± 0.03 
±0.03 
±2.03 
± 0.03 
±0.03 
±0.46 
±0.02 
±0.07 
±0.06 
± 0.04 
± 0.07 
± 0.03 
±0.46 
±0.03 
±0.05 
±0.03 
±0.01 
±0.03 



1.9: 
6.3: 
3.5: 
2.7: 
6.5: 
5.4: 
17.2: 
1.8: 

667.3 : 

1031.0: 

0.29±0.02 
28.5 ± 0.6 
182.5 : 

53.1 : 
75.3 : 
25.9: 
7.7: 



: 0.4 
:0.4 
: 0.4 
:0.4 
:0.4 
:0.4 
:0.5 
:0.5 
:2.9 
: 3.2 



13.2 
691.5 
7.8 
279.1 
10.9 
9.2 
6.2 
6.2 
640.6 
7.8 
1.5 
58.8 



: 1.5 
: 1.0 

:0.7 
:0.6 
:0.5 
:0.5 
:0.5 
:2.3 
:0.5 
: 1.4 
:0.4 
:0.4 
:0.5 
:0.5 
:2.4 
:0.5 
:0.2 
:0.8 



2.9 ±0.5 

19.8 ±0.5 

29.9 ± 0.6 
10.3 ± 0.5 
19.3 ±0.4 

9.4 ±0.5 
63.1 ±0.7 

6.8 ±0.5 
10.5 ± 0.4 
5.1 ±0.4 

7.5 ± 0.4 
5.4 ±0.4 



4.8 ± 0.3 
3.3 ±0.2 

4.6 ± 0.2 

3.5 ±0.2 
3.0 ± 0.3 
3.3 ±0.2 

20.4 ± 0.4 

1.2 ±0.2 
815.8 ± 3.2 

1791.0 ±5.4 
2.4±0.2 

15.0 ±0.3 
144.1 ± 1.4 
100.6 ± 1.5 

93.1 ±0.9 

25.5 ± 0.4 

7.8 ± 0.4 

4.3 ±0.2 
14.5 ± 0.5 

729.8 ± 2.3 

6.3 ± 0.4 
310.3 ± 1.5 

5.7 ± 0.2 

4.6 ±0.2 

2.8 ±0.2 

5.4 ± 0.3 
657.5 ± 2.4 

4.2 ± 0.2 

4.2 ± 0.2 
31.0 ±0.3 

1.3 ±0.2 
10.7 ± 0.3 

35.0 ± 0.3 
5.3 ±0.3 

11.1 ±0.2 

4.9 ± 0.2 
35.0 ± 0.4 

3.5 ±0.2 

5.9 ± 0.2 
3.0 ±0.2 

7.3 ±0.2 

3.4 ± 0.2 



91.6 ±0.2 
25.1 ±0.1 



1063.0 ± 1.0 
25604 ± 640 



4500 ± 2000 

602.9 ± 0.4 
54.4 ± 0.2 
101.7 ±0.2 



11272 ± 168 
1439.0 ± 0.6 



120.8 ± 0.3 



1.7 ±0.1 



64.3 ± 0.2 



E/1 

C/2AVB54 
D/3/MM 3 



HD 37903 
B/5AVB55 
1/7 

G/9AVB56 



12/W217/WB58 
14AV219 
11AV211 
F/15AV220 



11.4 ±0.2 
458.7 ± 0.3 H/16 



2.4 ± 0.2 



Kiso-A 



II? 
S 
I 
II 

II/III 
II/III 
II 
S 

I/II 
I 

II? 

s 

I/II 

I/II 

I/II 

II 

II 

s 

n 

I 

II 

I 

s 

s 

s 

n 

II 

s 

s 

s 

s 

s 

n 

s 

s 

s 

s 

s 

s 

s 

I/II 
s 



o 
o 



O 

2; 
o 
n 



" In column 1 2, integers coiTespond to NIR sources detected by Depoy et al. 1 1990), alphabet B-F are for sources identified bv lSellgrer] 119831) . WB stands for sources detected bv lWeaver & BabcocU )2004l) . W are 
sources detected bv lWitt. Schild & KraimanI U984.) and B33 names are from iBowler et al.j t200^) 



